Astronomy & Astrophysics manuscript no. 07gr sep21 djh print 


© ESO 2009 


September 21, 2009 





Extensive optical and near-infrared observations of tlie nearby, 
narrow-lined type Ic SN 2007gr: days 5 to 415 

D. J. Hunter^ S. Valenti\ R. Kotak', W. P. S. Meikle^, S. Taubenberger^ A. Pastorello\ S. Benetti^, V. Stanishev'^, S. 
J. Smartt', C. Trundle', A. A. Arkharov^-'', M. Bufano"^, E. Cappellaro', E. Di Carlo'^, M. Dolci', N. Elias-Rosa^ S. 
Frandsen^, J. U. Fynbo'", U. Hopp'' '^, V. M. Larionov^'^, P. Laursen'", P. Mazzali^'"^''^, H. Navasardyan'*, C. Ries'', 

A. Riffeser", L. Rizzi'"*, D. Y. Tsvetkov'^ M. Turatto'^, S. Wilke" 

' Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, BT7 INN, United Kingdom 
- Astrophysics Group, Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 lAZ, United Kingdom 
' Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str 1, 85741 Garching bei Miinchen, Germany 

INAF - Osservatorio Astronomico, vicolo dell'Osservatorio 5, 35122, Padova, Italy 
^ Pulkovo Central Astronomical Observatory, Pulkovskoe shosse 65, 196140, St. Petersberg, Russia 

* Astronomical Institute of St Petersburg State University, Universitetskij Prospect 28, Petrodvorets, 198504 St Petersburg, Russia 
' INAF Osservatorio Astronomico di CoUurania, via M. Maggini, 64100 Teramo, Italy 

Spitzer Science Center, California Institute of Technology, 1200 E.California Blvd., Pasadena, CA 91 125, USA. 
' Institut for Fysik og Astronomi, Arhus Universitet, Ny Munkegade, Bygn. 1520, 8000 Arhus C 

Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen, Denmark 
" Universitats-Sternwarte Miinchen, Scheinenstr I 81679 Miinchen, Germany 

Max-Planck-Institut fiir Extraterrestrische Physik, Giessenbachstr, 85748 Garching bei Miinchen, Germany 

Scuola Normale Superiore, Piazza Cavalieri, 7, 56126 Pisa, Italy 

Joint Astronomy Centre, 660 N. Aohoku Place, Hilo, HI 96720 
" Sternberg Astronomical Institute, University Ave. 13, 119992 Moscow, Russia 



Received / Accepted 



ABSTRACT 



We present photometric and spectroscopic observations at optical and near-infrared wavelengths of the nearby type Ic super- 
nova 2007gr These represent the most extensive data-set to date of any supernova of this sub-type, with frequent coverage from 
shortly after discovery to more than one year post-explosion. We deduce a rise time to 6-band maximum of 11.5 ± 2.7 d. We find 
a peak S-band magnitude of M^ = -16.8, and light curves which are remarkably similar to the so-called 'hypernova' SN 2002ap. 
In contrast, the spectra of SNe 2007gr and 2002ap show marked differences, not least in their respective expansion velocities. We 
attribute these differences primarily to the density profiles of their progenitor stars at the time of explosion i.e. a more compact star for 
SN 2007gr compai'ed to SN 2002ap. From the quasi-bolometric light curve of SN 2007gr, we estimate that 0.076 ± 0.010 Mq of '*Ni 
was produced in the explosion. Our near-infrared (IR) spectra clearly show the onset and disappearance of the first overtone of carbon 
monoxide (CO) between ~ 70 to 175 d relative to S-band maximum. The detection of the CO molecule implies that ionised He was 
not microscopically mixed within the carbon/oxygen layers. From the optical spectra, near-IR light curves, and colour evolution, we 
find no evidence for dust condensation in the ejecta out to about -1-400 d. Given the combination of unprecedented temporal coverage, 
and high signal-to-noise data, we suggest that SN 2007gr could be used as a template object for supernovae of this sub-class. 

Key words. Supernovae; general - Supernovae: individual (SN 2007gr, SN 2002ap) 



1. Introduction 

One of the main objectives of supernovae (SNe) research is to 
understand the relation between the physics of the explosion and 
the nature of the progenitor star and its immediate environment. 

Type I SNe, i.e. those lacking spectroscopic signatures of hy- 
drogen, are divided into the la, lb, and Ic categories. Of these, 
type la SNe constitute a physically distinct class and are be- 
lieved to arise from the thermonuclear explosion of an accret- 
ing carbon-oxygen white dwarf. Instead, both the type lb (He- 
rich) and the type Ic (He-poor) SNe are believed to originate 
from the core-collapse of stars more massive than ~ SMq. SNe 
of type lb and type Ic are quite heterogeneous in their spectral 
properties and energetics, which range from the relatively low 
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kineti c energy (~10^' erg) events like SN 19941 ( iNomoto et al.l 
11994 . to the high kinetic energy (~10^~ erg) broad-lined 
events sometime s dubbed 'hypernovae' e.g. the broad-lined type 
Ic SNe 1 998bw jGalama et al.lll998l: llwainoto et alJ [T998l) and 
2003dh jHiorth et al.l 120031: iMazzali et all l2003b . These latter 
SNe have be en associated with long-d uration gamma-ray bursts 
(GRBs) (see lWooslev & Bloomll2006l for a recent review). 

At present, there is no consensus regarding the differences 
between the progenitors that produce type lb and Ic SNe. 
However, it is generally accepted that Wolf-Rayet stars are 
the most promising candidates as these have shed their outer 
layer of H, as well as va rying amounts of their He layer 
dWheeler & Levreaulllll985l) . The H and He layers are primar- 
ily shed through strong stellar winds, resulting in WC and WO 
stars (i.e. carbon and oxygen dominated, respectively). The same 
mass-loss mechanism may likewise only remove the H layers so 
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that the potential progenitors of type lb SNe i.e. WN stars (ni- 
trogen dominating over carbon) exhibit the products of core-H 
burning. Alternativel y, the outermost layers ma y be stripped off 
by a companion star (iPodsiadlowski et aTll 1992b with the degree 
of envelope- stripping depending on the configuration of the bi- 
nary system dPols & Nomotolll997l) . 

The distinction between type lb and type Ic SNe is histori- 
cal, and challenged by the idea that He may be present in both 
types. The detection of spectral features due to He in type lb 
SNe might simply imply that the element is high in abundance 
and/or a sufficient quantity of radioactive ni ckel is mixed into the 



He layer, providing an excitation source ( 


Wheeler et al.ll 19871 


IShiaevama et al.ll 19901 iHachisu et al.lll991 


). The converse may 



be true for type Ic SNe expl aining the non-detection of optical 
He I l ines in type Ic spectra dNomoto et alJll990l : iHachisu et al.l 
1199 lb . Most important, no type Ic SN has yet displayed the iso- 
lated near-IR He line at 2.058 jjm to a strength comparable to 
those detected in type lb SNe. 

Consequently, it is possible that a continuous degree of 
He abundances and/or excitation exists from type lb to Ic 
SNe, making the distribution smooth rather than bimodal. This 
seems consistent with the identification of intermediate cases 
such as SN 1999ex which was characterised by weak op- 
tical He I Une s, and strong He I AA 10830, 20581 lines 
in the near-IR (iHamuv et al 1 12002|) . Additionally, the transi- 
tional nature of S Ne such as SN 2005bf (Polatelli et al. 2QQQ 
and SN 20080 (ISoderberg et al.1 120081: iMazzaU et al.1 l2008l 
iMalesani et al.ll2009l) . both of which underwent a metamorpho- 
sis from a type Ic at early times to a type lb at later epochs, may 
also provide an insight into the controversy of He. The current 
challenge lies in linking the observed variations of type Ib/c SNe 
to the physical properties of their progenitor systems. 

This paper presents the results of an intensive observational 
follow-up campaign of the type Ic SN 2007gr from shortly af- 
ter explosion to more than a year later The layout of this paper 
is as follows: in section |2] we provide some basic information 
on SN 2007 gr and its host galaxy, with the data acquisition and 
reduction procedures described in section 3. This is followed in 
section 4 by a description of the optical and near-IR photometric 
behaviour, including the colour evolution and quasi-bolometric 
light curve. Section 5 is devoted to the spectroscopic evolution; 
we end with a summary in section 6. 

2. SN 2007gr 

SN 2007gr was discovered during the course of the Lick 
Obs ervatory Supernova S earch (LOSS) on 2007 August 15.51 
UT dMadison & Li.ll2007l) . The SN exploded about 24'.'8 W and 
15'.' 8 N of the nucleus of the nearby spiral galaxy NGC 1058. 
IChornock et"an d2007h classified SN 2007 gr to be of type Ib/c 
from a spectrum taken soon after explosion. The classification 
was later refined to type Ic as th e presence of He co uld not be 
confirmed in subsequent spectra dValenti et al.ll2008al) . 

Given that the object was undetected in a n image taken 
on 20 07 August 10.44 (unfiltered mag < 18.9: iMadison & Li.l 
l2007h . the explosion date can be confined to less than 5 
days before discovery. Here, we adopt 2007 August 13 (JD 
2,454,325.5 + 2.5) as the explosion epoch, and 2007 August 25 
(JD 2,454,337.0 + 1.0) as the time of B band maximum, used 
hereafter as the reference epoch, t = days (see section 1431 ). 
Although no source was detected at the location of the SN in 
HST pre-explosion images (Crockett et al. 2008), the SN was 
suggested to be a cluster member based on its close location to a 
bright source. The pre-explosion HST and ground based images 



marginally favour a young cluster with a turn-off mass of 28 
+ 4M0. 

The host galaxy of SN 2007gr, NGC 1058, belongs to a 
group of nearby galaxies of which NGC 925 is also a member. 
The di stance to NGC 925 has been derived bv lSilbermann et all 
d 19961) as 9.29 + 0.69 Mpc (/z = 29.84 + 0.16) using Cepheid 
variables: we adopt this value throughout the paper for the dis- 
tance to SN 2007gr. We note that an alternative distance esti- 
mate based on the Tully-Fisher relation has also b een reported 
(10.6 + 1.3 Mpc: iTerrv. Paturel & Ekholmll2002l) . The galaxy 
of SN 200 7gr also hosted two pr ior supernovae, SN 1969L 
(Av = 017: 'Burstein & Heiles|11984|) and SN 1961V {Ay - 2.2- 
4. 1 : iFilippenko et al. . 1995ai) with th e latter possibly bein g the 
outburst of a luminous blue variable dGoodrich et al.l T989I) . 

We adopt a ga lactic extinction of E(B - y) = 0.062 
dSchlegel et alj [19981) . and a host galaxy extinction of E{B - 
y) = 0.03 which was estimated using the equivalent widths 
of Na I D absorptio n lines in a medium-resolution spectrum 
dValenti et alj|2008al) . We assume a total uncertainty of 20% for 
our estimate of the total extinction E(B - V),,,, =0.092 + 0.018 
mag. 

3. Data acquisition and reduction 

The proximity of SN 2007gr lent itself to an extensive observa- 
tional campaign starting 6 days before maximum light and ex- 
tending through to H-404 days. The data acquisition and reduc- 
tion procedures are detailed below. 

3.1. Photometry 

UBVRI photometry of SN 2007gr was obtained from 2007 Aug 
18 to 2008 Sep 30, spanning approximately -6 to +404 days 
with respect to B-band maximum. 

The basic reduction of the data i.e., trimming, bias and over- 
scan correction, and flat-fielding was performed using standard 
procedures within the IRAF environment. Images taken on the 
same night with the same filter were combined to improve the 
signal-to-noise ratio. 

The instrumental magnitudes of SN 2007gr were obtained 
via the PSF fitting technique as implemented in the SN00PM3 
software package. This technique was preferred over aperture 
photometry because of the complex location of SN 2007 gr, close 
to two bright sources (see inset Fig.[T]). 

Only for the observations taken on 2008 Feb 10, did we sub- 
tract a template image of the host galaxy before performing our 
measurements on the SN. This was necessary given the poor see- 
ing conditions (> r.'8). For the template, we used archival im- 
ages obtained on 2005 Jan. 13 with the Isaac Newton Telescope. 

The galaxy subtraction was performed usin g a purpose- 
written script based on the ISIS package dAlardI 120001 : 
lAlard&Luptonl 119981) and implemented into the SNOOPY 
package. Template subtraction was not applied to all the data 
as the only template images available, with similar filters to the 
Landolt system (Johnson- Cousins), were those from the Isaac 
Newton Telescope described above. The SN photometry was 
performed relative to a local sequence of 10 stars in the field 
of NGC 1058, calibrated during several photometric nights by 
comparison with Landolt standard stars (Fig. [1). The UBVRI 
magnitudes of the local standards are reported in Table I A. ll 

' SNOOPY was originally devised bv lPatatI ( 119961) and implemented 
in IRAF by Cappellaro. The package is based on DAGPHOT and has 
been optimized for SNe. 
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Fig. 1. Section of V-band image taken on 2008 Jan 13 (~ 
+ 142.4d relative to B maximum) with the Nordic Optical 
Telescope (NOT). SN 2007gr and the local sequence of stars are 
marked. The field of view is 6'. 3 x6'.3. The inset shows the two 
bright stars that flank SN 2007 gr, both of which are approxi- 
mately 5" from the SN. 



In spite of the many different instrumental configurations 
used, the colour-corrected UBVR photometry of SN 2007gr ap- 
pears to be coherent. In contrast, the raw /-band light curve 
showed a significant scatter (up to 0.4 mag) due to the differ- 
ences in the / filter transmission curves between the various in- 
strum ents. We therefore applied an S-correction to the / band 
data jPignata et al.l l2008l) in order to standardise the photom- 
etry. This technique has previously b een applied t o a number 
of typ e la SNe (e.g. lStritzinper et al.l 2002; Pignat a et al.ll2004 
|2008|) . Application of the S-correction to the /- band data sub- 
stantially reduced the scatter in the light curve. The S-corrected 
light curves of SN 2007gr are displayed in Fig.|2]and the magni- 
tudes are reported in Table A. 2. 

Near-IR photometry in the JHK bands was obtained on 25 
nights over a period of more than one year (Table IA.3l l. The 
near-lR magnitudes were determined using the same PSF-fitting 
technique as described for the optical photometry. The SN mag- 
nitudes were calibrated with respect to field stars in the 2MASSEI 
database (see section 14.21 for light curves). The derived JHK 
magnitudes and their uncertainties are listed in Table IA.3I The 
first-order colour corrections were omitted when observations in 
one band only were obtained. 

3.2. Spectroscopy 

Optical spectra were obtained at 29 epochs, spanning -7.3 to 
+ 375.5 days with respect to B maximum. Table lA4l lists the 
dates, instrument configurations, as well as the resulting spectral 
resolution and wavelength coverage. 

The 2D spectral images were processed in a standard fashion 
using IRAF, including bias subtraction and flat-field corrections. 
The spectra were wavelength calibrated using arc lamps; the 

^ http://www.ipac.caltech.edu/2mass/index.html 



caUbration was checked against night-sky lines and if required 
small shifts (< few A) were applied to the SN spectra. The spec- 
tra were flux calibrated using spectrophotometric standard-star 
spectra taken with an identical instrument configuration. All ob- 
servations were obtained with a slit angle chosen to minimise the 
contribution of the two stars flanking SN 2007gr (see inset Fig. 
[U. The slit was therefore not necessarily aligned along the paral- 
lactic angle. However, the spectra were corrected for light losses 
by making use of the extensive multi-band photometry available. 

Our last optical spectrum was taken on 2008 Sep. 3, approx- 
imately 387 days after explosion, with the Gemini Multi-Object 
Spectrograph on the Gemini North Telescope (GMOS-N). The 
reduction procedure for this spectrum differed slightly compared 
to the other spectra as we processed the data with the GMOS 
pipeline implemented within IRAF. The spectrum revealed sig- 
nificant contamination from nearby bright stars at bluer wave- 
lengths. The spectra of these nearby sources were extracted with 
the SN and their flux contributions were subsequently subtracted 
from the SN spectrum using standard tasks within IRAF. The op- 
tical spectra of SN 2007gr are discussed in section |5T| 

SN 2007gr was first observed in the near-IR on 2007 Sep 9 
UT (+ 15.5 d), and subsequently ten more times. The bulk of our 
near-infrared spectroscopic data were acquired using the near- 
IR imager and spectrometer (NIRI) mounted at the Cassegrain 
focus of the GEMINI-North Telescope. Three epochs were ob- 
tained using the Near Infrared Camera Spectrometer (NICS) at 
the TNG. Details of these observations, including instrument 
configurations, are given in Table lA.5] 

The NIRI data were processed using the GEMINI/NIRI 
pipeline implemented within the IRAF environment, while the 
TNG data were processed using standard routines in IRAF. All 
of our near-IR spectra were acquired by nodding the point- 
source targets along the spectrograph slit in the standard ABBA 
pattern, facilitating the removal of background emission. For 
the TNG data, the sky subtraction was performed by subtract- 
ing the A frames from the B frames (and vice-versa), and one- 
dimensional spectra of SN 2007gr were then extracted from the 
sky-subtracted frames. In the case of the Gemini data, the AB 
frames were pairwise subtracted and divided by the normalised 
flat field. A sky frame was then derived for each individual expo- 
sure by averaging the nearest 6-10 dithered frames taken within 
2 minutes. At later epochs, the time period was set in such a way 
so that one sky frame was selected per object image obtained on 
2007 Dec 20 and 2008 Jan 09. 

The near-IR spectra were all wavelength-calibrated using arc 
lamps. We used spectra of standard stars, taken at similar air 
masses and on the same nights as the SN spectra, to remove tel- 
luric absorptions and to obtain a relative flux calibration. The 
final flux calibration of the spectra was achieved by applying 
scaling factors derived from the near-IR photometry (Table |A3l l. 
Our series of near-IR spectra are shown in section |53] 

The spectra presented are in the galaxy rest-frame assuming 
a redshift of the galaxy NGC 1058 i.e. z= 0.001728 B From 
our latest optical spectrum on 2008 Sep 03 (+ 375.5 d), we mea- 
sure a redshift from the Ha line due to NGC 1058 that is con- 
sistent with the value reported above. The spectr a are corrected 
for ext inction using the standard extinction law of ICardelli et al.l 
(119891) and E{B - V),o, = 0.092 mag (see section|0]l. 



3 The redshift of NGC 1058 is taken from NED (Nasa/Ipac 
Extragalactic Database). 
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Fig. 2. Complete UBVRI light curves of SN 2007gr. /-band data 
obtained from the Liverpool and Wendelstein Telescopes have 
been S-corrected ( i]3.1l ). All errors at phases < 200 d are smaller 
than the symbol sizes. The dashed line shows the rate of decay 
of ^^Co— >^^Fe. The arrow indicates the pre-explosion/?-band de- 
tection limit at the location of SN 2007 gr The inset shows the 
light curve evolution at early times with the dotted line approxi- 
mately linking the peak magnitudes in all bands. 



4. Photometry 

4. 1. Optical Light Curves 

The complete optical light curves of SN 2007gr, extending from 
-6 to 404 d after B maximum, are displayed in Figure [2] The 
dense photometric coverage of SN 2007gr allows us to estimate 
the epoch of maxima and the peak apparent magnitudes in all 
bands. This was done by fitting each light curve around max- 
imum with a low-order polynomial function. Table [T] lists the 
peak apparent and absolute magnitudes, as well as the epochs of 
maxima in all bands with respect to both Z?-band maximum and 
the date of explosion. All pho tometry has been cor rected for ex- 
tinction using the relations of ICardelli et al.l (Il989h and a colour 
excess E(B - V)tor = 0.092 mag (see Table[T]l- 

Z?-band maximum occurred on 2007 August 25 at 13.47 
+ 0.02 mag, with the V, R, and / bands peaking respectively 
about 1 .5, 4.0, and 4.0 d later. With a distance modulus ju = 29.84, 
we estimate the absolute B-band magnitude of SN 2007gr at 
maximum to be -16.75. This is the same as for SN 2002ap 
{Mb - -16.76), but significantly fainter than the GRB asso- 
ciated SN 1998bw (Mb = -18.44) and SN 2003jd (Mb = 
- 19.30). Q Since the peak brightness depends primarily on the 
mass of ^^Ni synthesized in the explosion, we may immediately 
surmise that SNe 2002ap and 2007 gr appear to have comparable 
quantities of ^^Ni. We will discuss this further in section |4!4l 

The rise time of SN 2007gr, hence likely the time inter- 
val between explosion and maximum is constrained by the R - 
band non-detection on 2007 August 10 dMadison & Li.ll2007[) . 
This restricts the rise time to B maximum to 1 1.5 + 2.7 d (i.e. 
8.8- 14.2 d), indicating that SN 2007gr has an intermediate rise 
time compared to ot her type Ic SNe. It is typical to those 
of SN 2002ap (8d; iFolev et al.l l2003l) . SN 1998bw (-15 d; 



References are given in Table[2] 



iMazzaU et"ani200lh . and SN 19941 which reac hed B maximum 
~12 days after explosion dlwamoto et al.|[T994l) . 

We estimated the decline rates for SN 2007gr by using linear 
fits to the light curves (see Table[T]i. The UBVRI light curves re- 
veal a rapid decline after maximum followed by a much slower 
decline, which is typical of SNe Ib/c. The decline rate in the pe- 
riod from 30 to 170d after the explosion (i.e. ~ 0.01-0.017 mag 
d"') is similar in all optical bands, but typically a bit faster in 
the redder bands. The decline in the earlier part (30- 100 d) is 
generally less steep than at later epochs (100- 170d) with the ex- 
ception of the R band. The R light curve initially declines faster 
than other optical bands and is followed by a phase of slower de- 
cline from 100 to 170 d. This is probably due to the development 
of forbidden [Ca II] + [O II] ~ 7300 A and [O I] AA 6300, 6364 
arising at ~80d, as in the case of typical SNe Ib/c. At epochs 
later than 60 days, the light curve slopes of SN 2007gr in all op- 
tical bands are found to be steeper than those expected when the 
energy source is ^''Co— >^*Fe decay and the trapping of y-rays 
is complete [0.98 mag (100 d)"'] (see Fig|2]i. This is probably 
indicative of significant y-ray escape. 

In Figure [3] we compare the UBVRI light curves of 
SN 2007gr with those of other well observed stripped core- 
collapse SNe Ic. These SNe represent an exhaustive sample of 
type Ic SNe with well-sampled light curves; luckily, these SNe 
also span a large range of properties, thereby probing different 
physical conditions. 

The light curves of the SNe have been scaled in time and 
magnitude to match the peaks of SN 2007gr in all filters. Figure 
|3] illustrates the diversity of SNe Ic in te rms of their light curve 
morph ology and decline rates. Although IClocchiatti & Wheeler! 
(Il997h suggest that the decline rates appear to come in slow and 
fast varieties, new events seem to support the idea that there is 
a continuous range between these two extremes. Fig. [3] shows 
that after maximum brightness, SN 19941 is the fastest decliner 
in all filters followed by SN 2003jd, while SN 2004aw is gen- 
erally the slowest. In the nebular phase, the decline rates (ytc) 
of the light curves are determined by the energy deposition of 
the y-rays produced in the ^^Ni^^''Co— >^''Fe radioactive decay 
chain. Thus the emitted luminosity includes the y-ray contribu- 
tions which result from the energy deposition in the ''^Co decay 
and from the annihilation of electron-positron pairs, in which the 
kinetic energy of the positrons is an important factor Assuming 
homologous expansion and the incomplete trapping of y-rays 
and positrons, the behaviour of the light curve depends on the 
mass ejected (M^j) and th e kinetic ene r gy (£fa ,i) according to 
the following relation (see IValenti et al.l (l2008bh and references 
therein for a simple derivation), 

jLC^M.jlEun"^ (1) 

Hence, fast decliners such as SN 19941 (y^c ~0.88) are typical 
of SNe with lower ratios of ejecta mass to kinetic energy than 
the slow decUners of SN 2004aw (yz.c~ 2.3) and SN 1998bw 
(jLc ~ 1-8) (see Tabled. Meanwhile, the U band of SN 2007gr 
fades at a rate similar to those of SNe 2004aw and 2002ap from 
B maximum to ~50d later, while the later evolution is slower 
than in SN 2002ap. For the B, V and R bands, the decline is very 
similar to SNe 2002ap and 1998bw, and faster than SN 2004aw. 

4.2. Near-IR Liglit Curves 

Our near-IR photometry for SN 2007gr is reported in Table |A. 3] 
Early observations in the /-band enable us to constrain the peak 
magnitude and epoch of J maximum. This is likely to have oc- 
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y 






J 


H 




Apparent mag. 


13.29 ±0.02 


13.47 ± 0.02 


12.91 ±0.01 


12.76 ±0.02 


12.61 ± 0.02 


12.40 ±0.18 


12.06 ±0.10 


11.66 ± 0.10 


at maximum 




















Absolute mag. 


-16.99 ±0.19 


-16.75 ±0.19 


-17.22 ±0.18 


-17.29 ±0.18 


-17.36 ±0.17 


-17.52 ±0.20 


-17.83 ± 0.20 


-18.21 ±0.20 


at maximum " 




















Extinction (^mag} 


n 44 + n 09 


A -Jy . A A« 


29 + 06 


21 + 04 


1 3 + 03 


0.08 ± 0.02 


0.05 ± 0.01 


03 + 01 


JD ot 


j4jj4.0 ± U.o 


J^JJ / .u ± 1 .u 


±1.1 


j4j41 .U ± 1 . / 


j4j41 .U ± I . J 


54341.4 ±2.3 


54342.4± 4.3 


j4j4Z.J ± J.J 


maximum 




















Epoch of 


-2.2 ± 1.3 


0± 1.0 


1.5 ± 1.5 


4.0 ± 2.0 


4.0 ± 1.8 


4.4 ± 2.5 


5.4 ±4.4 


5.3 ±3.6 


maximum (days} 




















Rising 


9.3 ± 2.8 


11.5 ± 2.7 


13.0 ± 2.9 


15.5 ±3.2 


15.5 ±3.1 


15.9 ±3.0 


16.9 ±5.1 


16.8 ± 4.4 


times (tiays^ 




















Decline rates ^ 




















(30-lOOd) ' 


0.0095 ± 0.0007 


0.0111 ± 0.0006 


0.0170 ±0.0005 


0.0194 ±0.0007 


0.0164 ± 0.0004 


0.0304 ± 0.0006 


0.0268 ± 0.0008 


0.0342 ± 0.0008 


(100-170d) ^ 


0.0169 ± 0.0011 


0.0153 ± 0.0003 


0.0189 ±0.0003 


0.0128 ± 0.0004 


0.0167 ± 0.0005 










(30-170d) ' 


0.0114 ±0.0006 


0.0130 ± 0.0004 


0.0171 ±0.0003 


0.0159 ± 0.0005 


0.0169 ± 0.0002 














A distance modulus 


;i = 29.84 ±0.16 mag 


and a colour excess E(B-V) = 0.092 mag 


were adopted. 






^ The uncertainty in 


the eooch of maximum is defined as the time neriod in which the masnitude is not more than 0.01 mae fainter than the Dcak magnitude lYoshii et 


alJ2003l). 


' Epoch with 


respect to the estimated ZJ-band maximum JD 2,454,337.0 ±1.0. ^ Epoch with res 


pect to the estimated explosion 


date JD 2,454,325.5 ± 2.5. 








Decline rates (mag d ' ) were estimated using linear least-square fittinj 










Table 2. The parameters of type Ic SNe 
















SN (Ic) '"Ni mass 


Ejecta mass 


Ei-,„ mfl,„„„ 


Mbiuo," Rise 


E(B-V)" 


f 


Redshift' 


Ref 






(Mq) 


(Mo) 


(lO^' ergs) 


times (d) 










20071 


II 0.076 ± 0.020 2.0-3.5 


1-4 13.47 


-16.75 11.5 ±2.5 0.092 


29.84 


0.001728 


1,2 




2004aw 0.30 ± 0.05 


3.5-8.0 


3.5-9.0 18.06 


-17.63 


0.370 


34.17 


0.015914 


3 




20O3jd 0.36 ± 0.04 


2.5-3.5 


5-10 15.75 


-19.30 -13 


0.144 


34.46 


0.018860 


4 




2002ap 0.09, 0.07 


2.5-5.0 


4-10 13.11 


-16.76 -8 


0.090 


29.50 


0.002192 


5, 6, 7, 8, 9 




1998bw 0.7, 0.4, 0.5 


10.9 


20-50 14.09 


-18.44 -15 


0.060 


32.28 


0.008670 


10, 11, 12, 13 




1997ef 0.13 


9.6 


17.5 17.45 


-16.18 


0.000 


33.63 


0.011688 


14 




19941 0.07 


0.88 


1 13.77 


-17.06 -12 


0.300 


29.60 


0.001544 


15, 16, 17 





References: ' I Valenti et ai] I2008al) - This paper ^^ iTaubenberger et all 120061) ' TValenti et aUtlOOSH) ^ lYoshii et al.l 120031) '^ iFolev et all 12003ft ^ ICal-Yam et alj|2002h iTomita et all 
12006) ^ Mazzali et al.' ('2002) Iwamoto et al.'n998') " IGalama et al] fT99l) " INakamura etaTlfSOOl " IPatat et aljt200ll) " IMazzaU et alHSOOOl) " INomoto et"anfT994l) " 
IRichmond et aUU996i) '^ .Sauer et aU^2006) ^' The peak absolute magnitudes were estimated using the apparent magnitudes at maxima. Total extinction Redshift of the host 

galaxy (taken from the Nasa/IPAC Extragalactic Database NED). 



curred on 2007 August 29, ~16 days after the explosion date and 
4.4 days after B maximum. 



In Figure m we compare the JHK lig htcurves of SN 2007gr 
with those of SN 2002ap dTomita et al.l 12006) . which has the 
best coverage of all SNe Ic at near-IR wavelengths. In the near- 
IR, the light curves of SN 2007gr are very similar to those of 
SN 2002ap in terms of light curve morphology. Approximate 
peak magnitudes have been derived using the additional cover- 
age of SN 2002ap as a template, particularly for the H and K 
bands (see Table |2]i where the rising branch is unconstrained. 
The estimated peaks of SN 2007gr suggest that the time delay 
between B and J maximum is comparable to SN 2002ap (~5 
days for SN 2007gr and 7 days for SN 2002ap). The estimated 
decline rates for the period from +30 to + 140 d after B maxi- 
mum in the J, H and K bands are 0.029, 0.024 and 0.026 mag 
d"', respectively, again similar to those of SN 2002ap. However, 
the steeper decline rates of SN 2002ap at later epochs (> + 1 20 d) 
indicate a greater y-ray leakage from the radioactive decay than 
for SN 2007gr 



4.3. Colour Evolution of SN 2007gr 

The colour curves of SN 2007gr were produced using the optical 
and near-IR photometry detailed above. In Figure |5] we present 
theU-B,B-V,V-R,V-J,V-KandH-K colour curves 
of SN 2007gr and comparisons with the same sample of type Ic 
SNe: 2004aw, 2003jd, 2002ap, 1998bw and 19941. All curves 
have been corrected for reddening using the relations given by 
ICardelU et al.l ( |1989[) [see Table|2]for values and references]. 

The colour evolution of all SNe presented become monoton- 
ically redder from B maximum to ~ + 20 d. After this time, the 
differences between the SNe become evident. 

The early U-B colour of SN 2007 gr (Fig.|5l top-left panel) 
shows a steep rise by ~ 1 .5 mag. This increasingly redder colour 
may be due to both the cooling of the continuum and the pro- 
gressive increase of the metal line blanketing in the U band. 
After the peak, the U-B colour of SN 2007gr gradually turns 
bluer, followed by a plateau for ~ 60 d. Intermediate-phase data 
for SN 2002ap are lacking, but at t = 150 days it appears to have 
reddened considerably by ~ 0.3 mag with respect to the colour at 
~ + 20 d. SN 2007gr also marginally reddens around this epoch. 
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Fig. 3. Optical light 



curves of SN 2007gr compared to those of the Type Ic SNe 2004 aw dTaube nberger et al.' '2006b, 



2003id (lValentietal.ll2008b[) . 2002ap (iPandev et all 120031: iFolev et al.l 120031; iTomita et alj |2006|). 1998bw (.Galama et al 1998 
McKenzie & Schaefej 119991) and 19941 (lYokoo et al.l 119941: iGvoon et al.l (l995l: iLee et alJ 119951; iTsvetkov & PavlvukI IT995 



Richmond et al.lll996h . All the light curves presented have been shifted in time and magnitude to match SN 2007gr at maximum. 
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Fig. 4. The near-IR hght curves of SNe 2007gr and 2002ap 
(iTomita et al.ll2006h . The filled symbols joined by the full lines 
represent the data for SN 2007 gr. The data for SN 2002ap (open 
symbols, dotted lines) have been shifted in magnitude by 0.7, 
0.7 and 0.5 mag to match SN 2007 gr at ~ 20 d in the JHK bands, 
respectively. 



The B - V colour indices of all SNe (Fig. |5] , top-middle 
panel) in this sample rise at an approximate linear rate for the 



first 3 weeks before the B band maximum. After the colour peak, 
they rapidly become bluer, except for SN 2004aw which remains 
the reddest at 1.2 mag (~ +50d) before showing any sign of 
a decline. In the case of SN 2004aw, it is difficult to ascertain 
whether this behaviour is intrinsic to the SN, or simply due to 
an underestimate of the reddening as has already been suggested 
by lBenetti et al.l (l2006l) . On the contrary, SN 19941 remains the 
bluest SN (B - y ~ 0.2-0.5 mag at ~ H-50d) throughout its 
evolution. The B - V colour of SN 2007gr is intermediate be- 
tween these two objects, with strong similarity to SNe 1998bw 
and 2002ap. Uncertainties in the reddening estimates and pos- 
sible contamination from nearby blue stars preclude a further 
discussion of the bluer colours of SN 19941 dTaubenberger et aTl 
2006). 

The early V - R colour evolution (< + 80 d) of SN 2007gr 
is comparable to the evolution of SN 2004aw with a fast rise 
to maximum at V - R ~ 0.5 mag after ~ H-20d. Following the 
peak, all SNe tend to become bluer for a short period until ~ 
+ 70 d. This can be explained in part by cooling due to ejecta ex- 
pansion, with the peak of the spectral energy distribution shift- 
ing towards longer wavelengths. After this period, SN 2007gr 
gradually reddens, following the evolution of SN 2002ap which 
reaches a V - R colour of ~ 1.3 mag at ~ H- 300 d. The onset of 
this reddening is consistent with the development of the forbid- 
den [O I] AA 6300, 6364 emission and the blend of [Ca II] AA 
7291, 7323 and [O II] AA 7320, 7330 (see Fig.©. 

In Figure jSj the V - J and V - K colour curves of all SNe 
become monotonically redder, reaching a peak at ~ H- 30 d. This 
is then followed by a rapid decline to the blue up to H- 150d be- 
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fore undergoing a further transition to a redder colour In the 
near-IR too, the colour evolution of SN 2007gr best resembles 
SN 2002ap. The V - 7 colours of SNe 2007gr and 2002ap red- 
den to ~ 1 .5 mag at ~ H- 400 d. 

Signatures of any dust that might have formed in the ejecta 
of SN 2007gr should be evident in the near-IR colour evolution 
at late times. Furthermore, the H - K colour evolution tracks the 
evolution of the emission due to the first overtone of the carbon- 
monoxide molecule at > 2.3 fim (see Fig. |5] § 15.5b . which may 
be a precursor to dust condensation in the ejecta. Indeed, the 
H - K colour curve of SN 2007gr tallies with the appearance 
and disappearance of the CO-feature between ~ H-70 to H-200 
days relative to B maximum. Such an obvious signature of CO 
is not apparent in the H - K colour curves of any of the other Ic 
SNe shown in Fig. [5] with the possible exception of SN 1998bw 
which shows a mild upturn at about + 60 d, although the near- 
IR coverage for all objects other than SN 2007gr is admittedly 
lacking at epochs of interest. For SN 2007gr, given the rela- 
tively blue H - K colour {Eh-k is 1.2 mag) at about + 400 d, 

we conclude that dust formation is unlikely to have occured in 
SN 2007 gr by this epoch. 

4.4. The Bolometric (uvoir) Light Curves and Physical 
Parameters 

We constructed an Optical near-InfraRed (UVOIR) bolometric 
light curve by integrating the UBVRIJHK broad band fluxes. 
To derive the integrated SN luminosity, the magnitudes were 
corrected for reddening, converted to a flux density at the ef- 
fective wavelength, and finally integrated using Simpson's rule. 
The integrated flux was then converted to a luminosity using the 
distance modulus as assumed in § |2]i.e. yu = 29.84 + 0. 16. 

The bolometric light curve of SN 2007gr is shown in Fig. 
|6] For comparison, we have also incl uded the bolometric light 
curves of the SNe Ic sample. Following lValenti et al.l(l2008al) . we 
applied a toy-model to the bolometric light curve of SN 2007gr 
in order to constrain the physical parameters of the explosion 
(see Table |2]i. These include the ejected mass (Mej), the nickel 
mass (Msfi^i), and the kinetic energy (Eki,,). The procedure is 
based on a two-component analytical model to account for the 
photospheric and nebular phases of expansion. During the pho- 
tospheric phase (t < 30 d past explosion) homologous expansion 
of the ejecta, spherical symmetry, and the loc ation of the radioac- 
tive ^^Ni exclusively in the core is assumed (lArnett et al. 1982). 
At late times (t > 60 d past explosion), the model includes the en- 
ergy contri bution from the ^^Ni— >^^Co^^ ^Fe decay, following 
the wo rk of [Sutherland & Wheeled ( 1 1 9841) and ICappellaro et al.l 
( Il997h . In the nebular phase the ejecta are optically thin, so that 
the incompl ete trapping of y-rays and p ositrons has to be ac- 
counted for (IClocchiatti & WheelerllI997l) . 

The modelling of the bolometric light curve suggests that 
the ejected mass M^j of SN 2007gr (2-3.5M0) is similar to 
SN 2002ap but the kinetic energy /sa,,, is considerably lower 
(by a factor of ~3). However, we do note that the param- 
eters for SN 2002ap (see Table |2|i have been derived via a 
more involved m odelling procedure (e.g. iMazzali et al.l l200"2t 
iFolev et al.ll2003l) . 

A comparison of the early quasi-bolometric light curves 
shows that SN 2007gr appears to be most similar to SN 2002ap 
in terms of peak luminosity and width. As alluded to previously, 
the similarity in peak luminosity (MboI = -16.9) implies that 
these SNe produced similar quantities of ^''Ni (~ 0.07-0.1 Mq). 
We estimate that the quantity of ^^Ni ejected in SN 2007 gr is 
0.076 + 0.020, similar to SNe 2002ap and 19941. Interestingly, 



44 - 
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Fig. 6. The uvoir bolometric light curve of SN 2007gr and other 
type Ic SNe. The inset shows the differences in the widths of the 
uvoir lightcurves of this sample of type Ic SNe. 



SN 2007gr also resembles SN 2003jd in the shape of the early 
bolometric light curve, although at a much lower luminosity, 
while the rise to peak luminosity is comparable to SN 19941. 
The rise time of stripped envelope SNe depends both on the ra- 
dius of the pr ogenitor prio r to explosion and on the quantity of 
ejected mass (lArnettlll996h . For reference, we list the rise times 
of a sample of Ib/c SNe in Table |2] collected from the litera- 
ture, but note that for all objects other than SN 2007 gr, these are 
observationally poorly constrained, and the reality of any differ- 
ences in the rise times cannot at this point be directly attributed 
to physical differences between the progenitor stars. 

SNe 2007gr and 2002ap have similar peak widths which 
are evidently broader than that of SN 19941 (see Figure |6] in- 
set). Such broad peaks are usually attributed to more massive 
SNe Ic progenito rs (e.g. lArnett et al.lll982l:[Mazzali et al.ll2001l; 
iFolev et al.lfeool . Since the peak widths indicate the rate of en- 
ergy diffusion out of the core, this suggests that the progenitors 
of SNe 2007 gr and 2002ap had relatively large masses. We also 
expect the mass ejected in these two events to be larger compared 
to SN 19941 but less than for SNe 1998bw and 2004aw. 

The bolometric light curve of SN 2007gr indicates a late- 
time (120 to 380 day) decline of y = 0.0142 + 0.0002 mag d"'. 
These decline rates are comparable to those of the broad-lined 
SNe 1998bw and 2002ap which are y = 0.0156 + 0.0001 mag d"' 
and y = 0.0162 + 0.0003 mag d"', respectively. The decline rates 
are faster than expected from the decay of radioactive ""^Co (i.e. 
0.0098 mag d '). If the decay of ^^Co is responsible for pow- 
ering the late-time light curve of SN 2007 gr, then this suggests 
that the steeper decline is due to a decreasing y-ray trapping in 
the SN envelope with time. 

In Table |2] we compare SN 2007 gr to other type Ic SNe 
which show a wide range of physical properties at the explosion 
epoch. We can see that SN 2007 gr is most similar to SNe 2002ap 
and 19941 in terms of the amount of ^*Ni synthesised (see ij |4.4t . 
It can also be noted that SN 2007gr appears to be similar to 
SN 19941 with respect to the ratio of kinetic energy to ejected 
mass. Thus SN 2007 gr may be considered as the link between 
broad-lined events like SN 2002ap and low energy events like 
SN 19941. We note however, that the parameters listed in Table 
|2]have been derived using a variety of methods. 
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Fig. 5. Optical and near-infrared colour evolution of SN 2007gr compared with other SNe Ic: SN 2002ap, SN 2004aw, SN 1998bw, 
SN 19941. Evolution of the U - B (top-left), B -V (top-middle), V -R (top-right), V -J (bottom-left), V - K (bottom-middle) 
and H - K (bottom-right) colours. The insets in two of the upper panels show a zoom-in o f the early colour evolution of SN 2007 gr 
compared to SN 2002ap. All curves have been corrected for reddening. The ICardeUi et al. (.1989.) law was used to estimate the 
extinction in the different bands. See text for references. 



5. Spectroscopy 

5.1. Optical spectra 

The spectroscopic evolution of SN 2007gr, from 7.3 days prior 
to maximum light to + 375.5 days after, is densely covered and 
is shown in Fig. [71 T he line identifications are based on those by 
IValentiet al.] ( l2008al) . 

The early spectra exhibit Ca II H&K, Mg II ^ 4481, and 
Na I D /l/l 5891, 5897 which are responsible for the P-Cygni 
absorption features observed at 3810-3850 A, 4390-4490 A, and 
5720-5760 A, respectively. The Si II 6355 A absorption can also 
be identified in the early spectra, however it fades around two 
we eks after maximum ligh t . 

IChornock etaTI (l20()7h classified SN 2007gr as a type Ib/c 
from a spectrum obtained the night following discovery while 
waiting for follow-up observations. Eventually, no strong ev- 
idence for 6678 and 7065 A He features which are typi- 
cally in type lb SNe ~ 7 days after maximum light were 
found. Consequently, SN 2007gr was classified as a type Ic. 
IValenti et al] (l2008al) studied the feature at 6450 A, which may 
be misidentified as 6678 A He I in early-time spectra, and sug- 
gest that a more plausible identification is C II 6580 A at a ve- 
locity of ~ 11, 000 kms"'. Indeed this feature disappears around 
maximum together with another feature at 7000 A attributable to 
C II 7235 A, thereby providing support for the Ic classification. 

Moreover, almost two weeks after B maximum, one of 
the strongest optical transitions of the C I ion (9095 A) ap- 
pears, reaching velocities of ~8000 kms"'. The spectrum at 



around three weeks display two additional C I absorptions (9094, 
9406 A) due to the extended coverage at longer wavelengths 
(IValenti et al.l l2008al) (see section 15.21 and for further discus- 
sions). 

The first ten spectra, extending from one week before B max- 
imum to one week after, are characterised by an initially blue 
continuum which reddens significantly with time. Following 
maximum light, the Ca II near-IR triplet and the Na I D lines 
become very prominent in the spectra. Other absorption fea- 
tures include Fe II (4924, 5018, 5169 A) and Sc II (5552A). 
Approximately one week after maximum, the absorption troughs 
centred at 4200 A and 4800 A show the W- shaped profile that 
is characteristic of many type Ic SNe at a comparable epoch. 
By this time, the flux in the bluer bands is significantly reduced 
which is consistent with the redder U-B and B-V colours (Fig. 

13. 

The spectral evolution of SN 2007gr is relatively fast with 
the disappearance of Si II 6355 A already two weeks after 
maximum. The relatively early appearance of [Ca II] + [O II] 
~ 7300 A and [O I] AA 6300, 6364 at about H- 55 d represents the 
onset of the nebular phase. The absorption features in the spec- 
troscopic evolution of SN 2007gr tend to drift towards the red 
with time, as expected. This is a result of the formation of lines 
occuring in slower, deeper layers of the ejected material. 
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Fig. 7. Evolution of the optical spectra of SN 2007gr. The spectra have been corrected for reddening, E(B - y),of = 0.092, and for 
the host galaxy redshift of NGC 1058, z = 0.001728. The spectra have also been magnified and displaced vertically for clarity by the 
numbers shown in brackets. Identifications are shown in the spectra of SN 2007gr at -6.5 and 158.4 days relative to B maximum. 
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The late-time spectra are dominated by the emission of for- 
bidden Unes such as Mg I] 4571 A, [O I] AA 6300, 6364, [Ca II] 
AA 7291, 7323, and the permitted transitions of Na I D, 1 7774, 
and the possible blend of C I (8727 A) and the Ca II near-IR 
triplet. 

The optical spectra do not reveal any clear indication for the 
presence of He I lines. If He I 5876 A was present it would al- 
most certainly be contaminated by Na I D. Other lines of He I in 
this region (e.g. 4471, 7065 A) would also be blended with other 
strong features at these wavelengths. Since the optical spectra 
alone are insufficient to detect small amounts of He I, in sec- 
tion 15.31 we examine the near-IR spectra of SN 2007 gr. The 
He I 1 0830 A Une is expected to be the strongest He I tran- 
sition jSwartz et al. 1993) and therefore s hould be highly sen- 
sitive to very sm all quantities of helium dWheeler at al.|[l993t 
iBaron et al.|[l996l) . 

5.2. Optical Spectra Comparison 

In Figures [8] (& |9l we present early and late-time spec- 
tral comparisons of SN 2007 g r with the SNe Ic sam- 
ple: SN 199 8bw (ISollerman et all 12000; Mazzali et al. 2001), 
SN 20 03id dValenti et al.| |2008bl). SN 20 02ap dGal-Yam et al] 
I2002t iKinugasa et al.l l2002 t iFolev et all l2003h. SN 2004aw 
dTaubenberger et al I l2006l) and SN 19941 (iFilippenko et al.1 
Il995bl : [Sauer et ai1l2006h . All of the spectra have been ordered 
according to their explosion energies (top spectra resulting from 
the most energetic) and corrected for reddening and redshift us- 
ing the values listed in Table |2] We have previously emphasised 
the remarkable similarity in light curve morphology and evo- 
lution between SNe 2007 gr and 2002ap. We will see that this 
similarity does not extend to the spectral line profiles and ejecta 
expansion velocities. 

The near-maximum spectrum of SN 2007gr resembles that 
of SN 19941 on day H-1 in terms of the presence of all major fea- 
tures. The main differences between these events at this epoch 
appears to be the deeper absorption of O I 7774 A relative to Na 
1 5896 A and the narrower spectral widths in SN 2007gr. In con- 
trast, the spectra of SNe 1998bw and 2002ap exhibit very broad 
features as a result of line blending from high-velocity ejected 
material. The narrow lines of SNe 2007gr and 19941 allow for 
a more secure line-identification and a better separation of spec- 
tral features. This is especially true for the C I absorptions (8335, 
9094, 9406 A) which are evident in SN 2007gr at H-20.6d and 
are also detected in SN 19941, although at weaker intensities. 
For SNe 2003jd and 2004aw, it is possible that the individual C I 
lines may be disguised as a single absorption blend at ~ 9000 A. 

SNe 2007 gr and 19941 also have better developed P Cygni 
profiles for the Na I D and Ca II near-IR triplet lines at only two 
weeks past maximum as a consequence of their rapid evolution. 
The comparatively slowly evolving SNe 1998bw, 2004aw, and 
2003jd, show no sign of these features until ~2 months. 

Four months after B maximum, the spectra of SN 19941 are 
the most evolved and show only forbidden line emissions. The 
absence of the photospheric O I 7774 A emission line also indi- 
cates that the evolution from the photospheric to nebular phase is 
faster than in SNe 1998bw, 2002ap and 2007gr. Even at this neb- 
ular phase, the spectral lines of SN 2007gr remain narrower than 
the other type Ic SNe presented, indicating the slower expansion 
velocities of the ejecta. The nebular phase spectra all show [O I] 
AA 6300, 6364 as one of their strongest features, complemented 
with lines of [Ca II] AA 7291, 7323 / [O II] AA 7300, 7330, and 
the near-IR Ca II triplet. The largely isolated nature of the [O I] 



doublet is found to be narrower in SN 2007gr compared to all 
the SNe Ic presented including SN 1 9941 by ~ 1 000 km s" ' . 

The nebular spectra of SN 2002ap, 1998bw and 2007 gr have 
strong [O I] AA 6300, 6364 and Mg I] A 4571 emissions. The 
ratio of Mg I] to [O I] may provide some information on the de- 
gree of stripping of the outermost layers experienced by the pro- 
genitor star, with a higher Mg I]/[0 I] ratio implying a greater 
degree of strippin g as more of the O-Ne-Mg layer is exposed 
(iFolev et al.ll2003l) . Additionally, a higher Mg I]/[0 I] ratio at 
later epochs may simply imply that less of the outer C/O layer 
is observed due to the transparency of the ejecta. We note that 
the [O I] and Mg I] lines may well be formed via different 
mechanisms, with [O I] being collisionally excited while Mg 
I] arising predominantly du e to recombination as suggested by 
iKozma & Franssoni (Il998l) . This could influence the line ratios 
that we measure in different SNe. Detailed modelling of the ex- 
citation mechanisms for these lines is beyond the scope of this 
work, so we merely flag this caveat, and use the line ratios as a 
proxy for the composition. 

In Figure [TO] we show the evolution of the Mg I]/[0 I] ra- 
tio for a sample of CCSNe during the nebular period H- 100 
to H-400d. The Mg I]/[0 I] ratio was calculated by first fit- 
ting a cubic spline to a local continuum. The continuum was 
then subtracted so that the integrated flux in the emission lines 
could be estimated. All of the SNe in Figure [TO] that have more 
than one epoch of measurement available, show a trend of in- 
creasing Mg I]/[0 I] with time, with the marked exceptions of 
SNe 1998bw and 2006aj. Interestingly, these two SNe are the 
only SNe in this sample that were associated with long-duration 
GRBs. For the SNe for which only a single measurement was 
possible (SNe 1985F, 1997ef, 2000ew, 2003jd), the Mg I]/[0 I] 
ratio was found to lie within the bounds spanned by other SNe 
of this type, and are included in Fig.[TO]for completeness. 

The increasing Mg I]/[0 I] trend observed for the majority 
of the SNe sample may simply be du e to an abundance effect, 
as hypothesised by iFolev et akl (|2003|) . In symmetric or quasi- 
symmetric explosions, the onion-structure of the progenitor star 
is almost preserved. This is also true in non-symmetric explo- 
sions viewed equatorially. In these cases as time goes by we ob- 
serve deeper and deeper into the ejecta reaching the O-Ne-Mg 
shell, and so the Mg I]/[0 I] is expected to increase. 

At the beginning of the nebular phase {+ 120 to + 170 d), 
SN 2007 gr has lower Mg I]/[0 I] ratios than the other SNe Ic. 
The lower Mg I]/[0 I] ratio of SN 2007 gr at these early epochs 
may suggest that we are still observing the highly concentrated 
C/O layer of the outer ejecta. This may support the idea that 
SN 2007gr originated from a carbon-rich progenitor (see i^ l5.3b . 
Another possible reason for the lower Mg I]/[0 I] ratio may 
be due to the inefficient burning of the large quantities of car- 
bon in SN 2007gr, thus producing only small quantities of Mg. 
Conversely, the larger Mg I] to [O I] ratios of SN 2002ap may 
imply that the progenitor experienced greater stripping of the 
outer carbon layer, since at an epoch similar to SN 2007gr we 
are already probing deeper into the core of the p rogenitor, re- 
vealing the O-Ne-Mg layer (e.g. lFolev et 



: or tnepn 
an i2003h . 



The data also show that SNe 2007gr and 2002ap exhibit 
larger Mg I] to [O I ] ratio s than the type lb SNe 1985F 
('Fil ippenko & Sargent! Il986l) and 2007Y dStritzinger et al.1 
,20091) . In this figure, we have also included th e ratios of 
the type lib SNe 2 00 lig ([Silverman et al.1 l2009l) and 1993J 
(iFransson et al.l2005h . The Mg I]/[0 1] ratios of the type lib SNe 
are also consistently lower than those of the type Ic SNe 2007 gr, 
2002ap and 2004aw at the same epochs. 
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Fig. 8. Comparison of spectra of SNe close to B-band maximum (left) and two-three weeks post maximum light (right). All spectra 
are de-reddened and redshifted assuming the values reported in Table|2] Typical absorption lines are identified. The three solid lines 
in the right-hnd panel next to SN 2007gr indicate the positions of the C I lines. SN 19941 also shows the C I features, although these 
are not as strong as in SN 2007 gr, while the other SNe show a broad dip in this region; these variations are probably attributable to 
the differences in the respective ejecta velocities of the SNe shown here, with the C I features being smeared out in those with the 
highest velocities. 
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Rest Wavelength (A) 

Fig. 9. Comparison of spectra of SNe about one month (left) and four months (right) after B maximum light. The P-Cygni absorp- 
tions of Na I D (AA 5891, 5897), O I (A 7774) and the Ca II near-IR triplet are clearly marked in the left panel. Typical nebular 
emission lines are indicated in the right panel. 
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Fig. 10. The Mg I]/[0 I] ratio as a function of time for a sample of SNe. All measuremen ts were performed as d escribed in section 
15.21 with the exception of type lib SN 200 lig, where the measurements were t aken from Silverman et al. ( 20091). All other Ic SNe 
spectra are from references listed in Table |2] with the exception of SN 2 006aj ([Mazzali et alJ 20071 : iTaubenberger et aDl2009l) and 
2000ew (Valenti et al. in prep). T he spectra of the type lb SNe 2007Y jStritzinger et al.ll2009l) and 1985F dFilippenko & SargentI 
119861) . and the type lib SN 1993 J dFransson et al.ll2005i) have also been used. They were corrected for extinction before the mea- 
surements were made. 
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Fig. 11. Left : The profiles of [O 1] 6300, 6364 and Mg 1] 4571 at 158 d past B maximum. The profile of [O 1] is plotted in velocity 
space using a rest wavelength of 6315 A, as obtained via gaussian fitting from the modified profile of [O 1] in the normal and low- 
density limits, 1:1, 3:1. The Mg 1] profile includes the contribution of an artificial component with an intensity ratio of 0.33 relative 
to the real peak and the same velocity offset as found in the oxygen doublet. The peaks of the Mg 1] profiles have been scaled to 
match the peak of the oxygen doublet at + 158 d. Centre : Evolution of forbidden [O 1] 5577, Mg 1] 4571 and [O 1] 6300, 6364. 
Right : Evolution of the profiles of the permitted lines of O 1 7774 A and Mg 1 15025 A. The vertical line in the panels correspond 
to zero velocity i.e. the rest wavelengths of the line profiles. 
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The general trend from our sample supports the idea that the 
progenitors of type Ic SNe experience greater envelope stripping 
than those of type lb and lib. 

The decreasing Mg I]/[0 I] trend observed for SNe 1998bw 
and 2006aj may be explained by the rapid rotation of the GRB- 
progenitors which result in greater mixing (via meridional circu- 
lation), and probably induce some degree of asymmetry. A com- 
plete mixing of the onion structure is expected in non-symmetric 
explosions viewed along the polar axis (along the jet). However, 
it may be that because of the efficiency of this mixing mech- 
anism, the early Mg I]/[0 I] ratios will be close to the aver- 
age values observed in other Type Ic SNe (e.g. first point of 
SN 1998bw). As time progresses, the quantity of Mg will di- 
minish and since O is more abundant, its lines will dominate 
over Mg lines. In this simple view, the SNe which show the 
trend of SNe 1998bw and 2006aj could be jet-like explosions 
or explosions where the ejecta has been fully mixed. An alterna- 
tive interpretation may be that the ionisation equilibrium evolves 
with time for Mg and O (which have very different ioni sation po- 
tentia ls) simply because of decreasing temperature (cf. lMeverottI 
fT980l) . This may well be a contributing factor that might account 
for the observed differences. 

Recent studies of nebular spectra suggest that around 
half of all SNe Ib/c are strongly aspherical while the re- 
mainde r appear to have moderate signatur es of aspheric- 
ity (Maeda et al.1 120081: iModjaz et al. 2008; Taubenberger et al.l 
[200 9; Milisa vhevic et al.1120091) . As previously mentioned, [O I] 
AA 6300, 6364 and Mg I] /1 4571 are two of the most dominant 
features in the late-time spectra of stripped-envelope SNe. As a 
consequence of their relatively isolated nature and higher abun- 
dances to other nebular lines, an analysis o f the emission profile s 
permits one to study the ejecta geometry jMazzali et al.ll20()5l) . 
In particular, one might expect Mg and O to have a s i milar spa- 
tial distribution within the ejecta (e.g. Maeda et al. 2006 ) and 
thus the profiles of these emission lines should also be similar 

Unfortunately, a direct comparison of the Mg I] and [O I] 
lines is not possible since the [O I] doublet arises from two for- 
bidden transitions which share the same upper level (^Pi,2 - 'D2). 
The relative intensity of the A 6300 and A 6364 lines will depend 
on the O I density and may vary from 1:1 to the low-de nsity 
Umit of 3: 1 (ISpvromiliolll994l) . iTaubenberger et al.l (|2009|) sug- 
gested that the density in stripped envelope SNe should, after 
lOOd, have decreased eno ugh to fi x the line intensity ratio to 
3:1, but recently [Milisavljevic et al.l ( |2009|) suggested a ratio 1 : 1 
to explain the Ol AA 6300, 6364 profiles in 4 stripped envelope 
SNe. In order to investigate which line intensity ratio is more 
appropriate for SN 2007gr, we have displayed the profiles of [O 
I] AA 6300, 6364 and Mg I] A 4571 in Figure [TT] (left panel) 
at 158 d after maximum. Since the profiles of Mg I] {+ an arti- 
ficial second component with a line intensity ratio of 0.33 and 
the same velocity offset as found in the oxygen doublet) and 
[O I] are similar, this tells us that the ratio between O I 6300 
and 6364 is close to 3:1. We also performed a comparison be- 
tween [O I] AA6300, 6364 and [O I] i5577 {+ an artificial sec- 
ond component with ratio 0.33, 0.66, 1.00) without obtaining a 
good match. The [O I] emission profile contains one peak po- 
sitioned at a blueshift vel ocity of ~-1500 km s~'and the other 
skewed at ~-800 kms ITaubenberger et akl (l2009h found an 
average blueshift of ~ 1000 km s ' for the [O I] /f 6300 line cen- 
troids of a large sample of Type Ib/c SNe which is also consistent 
with SN 2007 gr. The authors suggest that the observed effect is 
probably related to the ejecta geometry and may be a signature 
of a one-sided explosion, or, more likely, residual opacity in the 
inner ejecta. 



We also present the evolution of the permitted and forbidden 
lines of O and Mg at phases > H- lOOd (centre and right panels 
of Figure fTTTi. The figures show that the permitted lines of O I /I 
7774 and Mg I A 15025 behave very differently than the forbid- 
den lines of [O I] A 5577 and Mg I] A 4571. The emission peaks 
of the permitted profiles are found to be almost centered at their 
rest wavelengths while those of the forbidden lines are strongly 
blueshifted, with the effect becoming progressively weaker with 
time. In addition, the [O I] profile at 5577 A appears to be even 
more strongly blue-shifted compared with other forbidden lines. 

The observed trend of blueshifted forbidden line centroids 
compared to the central permitted lines suggest that they appear 
in different regions of the ejecta since we know that the forbid- 
den lines emerge from a lower-density medium. The systematic 
effective blue-shift of the forbidden lines may well be due to the 
residual opacity of the intervening material such that the photons 
emitted on the rear side of the SN are absorbed along the line of 
sight; 

ITaubenberger et al.l (|2009|) and iMilisavhevic et al.l (|2009|) 
provide alternative explanations for the observed trend, includ- 
ing dust formation and contamination from other ions. However, 
in section |53] we claim that dust has not yet formed in the ejecta 
of SN 2007 gr Additionally, any dust that is present should actu- 
ally result in a line blueshift which increases with time, not the 
contrary. Contamination due to other emission lines may also be 
neglected due to the similarity of the [O I] and Mg I] profiles, 
while contamination on the [O 1 1 A 5577 cannot be e xcluded. 

For SN 2007 gr, we note that lTanaka et al.1 (|2008|) report non- 
zero polarization across the Ca near-IR triplet at ~21d post- 
maximum. They concluded that their observations could be in- 
terpreted either in the context of a bipolar explosion with an 
oblate photosphere viewed from a mildly off-axis line-of-sight, 
or a spherically symmetric photosphere with a clumpy Ca II dis- 
tribution. Out to about 400 d, the optical and near-IR spectra, 
show no evidence of double-peaked profiles as would be ex- 
pected for a significantly aspherical explosion. From our data, 
we conclude that the explosion of SN 2007gr could not have 
been more than mildly aspherical. 

5.3. Near-infrared spectra 

All near-infrared spectra of SN 2007gr are presented in Fig. [12] 
At early epochs, the shape of the NIR continuum is nicely ap- 
proximated by a black body continuum. The earliest spectra are 
almost featureless with few lines due to intermediate mass ele- 
ments such as Mg I, O I, Na I. After maximum light, lines due 
to the iron group elements begin to appear 

One of the most prominent features in the near-IR spectra 
is an absorption at ~ 1.04yum. The origin of this feature has 
been discussed widely b y several authors, espe cially in the con- 
text of type la SNe, with 'Wheel er at al.l (1 19981) favouring Mg II 
and Meikle et al. ( 1996) suggesting a blend of Mg II and He I 
1.083 fim; the latter is a line that is easily excited in a variety of 
astrophysical contexts. In fact, other elements such as C, O, Si, 
and Fe may also contribute to features in this region. 

If He I were indeed responsible for the absorption, it would 
only be detected if the lower level of He becomes highly popu- 
lated in the presence of high energy, non-thermal radiation. The 
source of this radiation would be the y-rays produced in the ra- 
dioactive decay of ^*'Ni to ^''Co. The lower level population may 
then absorb photospheric photons of ~1 micron yielding the P 
Cygni profile. 
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Fig. 12. The near-infrared spectral evolution of SN 2007gr The spectra ar e displa yed in the rest frame of NGC 1058. To increase 
the temporal coverage we include the spectrum obtained bv IValenti et al.l (12008 al) on day + 15. The + 376.6 d spectrum has been 
smoothed with a box size of 5 pixels. The spectra have been displaced vertically for clarity by the numbers shown in brackets. The 
rest wavelengths of possibly identified fines are also included. 
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lMazzaU&Lucvl(ll998l) were able to reproduce the 1 .05 mi- 
cron feature in the Type la SN 1994D using just 0.01 of he- 
lium, implying that only small quantities are needed to account 
for this feature. Thus the apparent lack of strong He lines in type 
Ic SN spectra may imply a low abundance of He I and/or insuf- 
ficient mixing of ^^Ni into the outermost layers of the ejecta. 

If He I were a major contributor to the 1.04/Ym feature in 
SN 2007 gr, this would imply an average expansion velocity of 
-13,000 km s"' at 15.5 d which later declines to -10000 km s"' 
by 1 14.5 d. In spite of the fact that SN 2007gr is a narrow-lined 
Ic, so that the identification of lines should be facilitated, the op- 
tical region is plagued by blends, especially in the regions where 
He I lines would lie e.g. the He I 3888 A line lies close to the 
Ca II H & K Unes and the He I 5875 A lies close to the Na I 
D. The possible presence of the He I line at 5875 A is suggested 
in Section 5.4 to account for the unusual velocity distribution of 
the Na I D feature. In the near-IR, one might expect to detect He 
I 2.058 lyL/m. However, we see no evidence of any feature at this 
wavelength in our entire near-IR spectral series. 

If He I were a major contributor to the 1.04/im feature in 
SN 2007gr, this would imply an average expansion velocity of 
-13,000km s^' at 15.5 d which later declines to -10000km s^' 
by 1 14.5 d. In spite of the fact that SN 2007gr is a narrow-lined 
Ic, so that the identification of lines should be facilitated, the op- 
tical region is plagued by blends, especially in the regions where 
He I lines would lie e.g. the He I 3888 A line lies close to the 
Ca II H & K Unes and the He I 5875 A lies close to the Na I 
D. The possible presence of the He I line at 5875 A is suggested 
in Section 5.4 to account for the unusual velocity distribution of 
the Na I D feature. In the near-IR, one might expect to detect He 
I 2.058 1/im. However, we see no evidence of any feature at this 
wavelength in our entire near-IR spectral series. 

Although He I 2.058 l/^m is not visible in SN 2007 gr, we 
cannot fully disregard the contribution of helium to the 1.04 
jum absorption. The longer decay time from the 2s3S level of 
He I 1 .0830 /im means that this line can behave much like a low- 
energy, ground state transition, unlike the faster (x 4.5x10^ s) de- 
cay time of He 1 2.058 l/^m, and so will be easier to excite. Thus 
the absence of the He 1 2.058 1 fim line does not necessarily mean 
that the 1.0830yL(m line does not contribute to the 1.04/Ym fea- 
ture. This cannot be quantified without detailed modelling which 
is beyond the scope of this paper Yet, we can conclude that if He 
is present in SN 2007gr, the amount must be small. 

IValenti et al.l ([2008a) suggested that the 1.04yum feature in 
SN 2007gr may be attributed primarily to the C I ion (rest wave- 
length at 1.0695/im). They found that earlier (optical) spectra 
showed evidence for C II features at very early epochs (- -8 d) 
which fade rapidly, to be replaced by lines due to C I. 

From an initial modelling effort, they obtained a fairly rea- 
sonable match to their optical-near-IR spectrum at day H- 15, also 
reproducing several other observed C I lines. However, even with 
a model enhanced in C I, the entire 1 .04 yum feature, and the blue 
wing in particular, could not be matched. They suggested that 
this could be due to the contribution of one or more other ions 
e.g. Mg II, Si II. Indeed lMiUM-d et al.1 ( [19991) also had difficulty 
in accounting for this infrared absorption feature in SN 19941 
with He I alone, suggesting that it could be a blend of He I and 
C I lines. 

Together with the 1.04//m feature, the nebular spectra of the 
J band reveal emission lines of Fe II, Si II, Mg I, Na I, O I 
and S I, while the H band is dominated by emission features 
produced by iron-group elements (Fe II). The K band spectra 
appear to contain emission features from Na I, Mg II and the 



first overtone of CO (see also Figure [T6]l. We will discuss this in 
The [Fe II] lines at Aresi= 1.257//m and Aresi= 1.644yL(m are 
particularly prominent in the J and H bands, respectively. These 
are complemented with unblended Fe II lines in the optical at 
expansion velocities - 4000-6000 km s"' . 

The broad feature positioned at -1.13 /im can be identified 
with the 1.1287;um (3p ^^P-3d ^^D °) O I transition dMeikle et all 
il989) . The Mg I (1.503yum) feature observed in the H band is 
matched by strong features of Mg II at 4481 A in the optical re- 
gion. The K band region of the NIR contains the 2.207 fim dou- 
blet of Na I. In the optical too, Na I dominates the region around 
5900 A from the early to very late phases. Features which appear 
near 1.94 and 1.98yum around 138 d can probably be identified 
with Ca I. 

Currently, few near-IR spectra of Ib/c SNe are found in 
the literature. This is particularly true at late phases. At earlier 
epochs, the situation is somewhat less dire and we compare our 
early time spectra of SN 2007 gr with those of SNe Ic 2004aw 
(Taubenberge ret al. 2006,), 1998b w (Patat et al. 2001), the SN lb 
1999ex dHamuv et al.l2002b . and previously unpublished spectra 
of SN 2002ap (Icfl The near-IR spectra of all SNe mentioned 
above are presented in Fig. [13] 

The broad-lined spectra of SNe 2002ap and 1998bw ap- 
pear almost featureless, with the exception of a weak absorption 
at -1.04jum. This behaviour probably reflects the severe line 
blending due to high ejecta velocities in these objects. In con- 
trast, the distinct narrow features of SN 2007gr are clearly vis- 
ible at similar epochs. The type lb SN 1999ex exhibits equally 
strong He I absorption features at both 1.04jum and -2.0 jum. 
This is supported by the p resence of He I Unes in the optical 
spectra (iHamuv et al.ll2002l) . Interestingly, all SNe shown in Fig. 
[T3]display a prominent feature near 1 .04 yum regardless of epoch, 
however with the exception of SN 1999ex none shows an absorp- 
tion near 2.0 yum at any time. The 1.04jum absorption feature in 
SN 2007gr appears to be stronger than in SN 1999ex, providing 
some support for the contribution of C I. 

5.4. Ejecta velocity 

The velocities of several ejecta lines were determined by fitting 
a Gaussian profile to the absorption component of the P-Cygni 
feature in order to measure the position of the blue-shifted min- 
imum. This method provides an estimate of the expansion ve- 
locities of the layers where the individual lines predominantly 
form and therefore the distribution of elements within the ejecta. 
Only the velocities of a few spectral features have been mea- 
sured since for many others line blending is a problem, even 
for the relatively low velocities of SN 2007gr The evolution of 
ejecta velocities is shown in Fig. [14] 

The Ca II velocities of both the H & K features as well as 
the near-IR triplet are significantly higher than those of Mg II 
(4481 A), Fe II (4924, 5018, 5170 A), Na I (5890-5896 A), O I 
(7774 A), C I (9087, 10695 and 11753 A). Assuming homolo- 
gous expansion, the higher velocity Ca II lines imply that they 
form further out in the ejecta than the O I, Na I, and Fe II lines 
which emerge from the inner ejecta layers. The formation of Ca 
II lines further out is predominantly because these lines are ex- 
tremely strong and therefore become optically thick very rapidly, 
even in the low density, outer layers of the ejecta. As expected, 
the absorption component of the P-Cygni minima gradually shift 



' These spectra were obtained at the UKIRT telescope and reduced 
using standard procedures. 
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Fig. 13. Comparison of the near-IR spectra of SN 2007gr to SNe 
2002ap (broad-lined Ic), 1998bw (broad-lined Ic), 1999ex(type 
lb), and 2004aw (type Ic). Note the presence of the prominent 
He I 1.083, and 2.0581 jum lines in the spectrum of SN 1999ex. 
See text for references. 

redwards with time revealing the deeper and more slow moving 
ejecta. 

Comparisons of velocity widths of the Na 1 A 5891, Si 11 
A 6355, O 1 A 7774 and the Ca 11 NIR triplet for SNe 2007gr, 
2004aw, 2002ap and 19941 are shown in Fig. [15] The expansion 
velocity Vexp of SNe Ic can be directly measured from the min- 
ima of spectral absorption lines, of which the Si II line (6355 A) 
is commonly used when present. Assuming that the blue mini- 
mum contains no contamination from other ion species, the de- 
rived Vexp values of SN 2007gr reveal comparatively low veloc- 
ities which are probably related to the narrow features observed 
in the spectra. Using the Si 11 velocity as an approximation, we 
estimate a photospheric velocity of ~ 6,700 kms"' at B maxi- 
mum, which decreases to only ~ 4,000 km s ' ten days later 

In contrast, the photospheric velocity of SN 2002ap declines 
rapidly and ranges from ~ 17,000 to 7,000 kms"' in a simi- 
lar time frame. The blue-shifts of the Si II line of these two 
objects are comparable at two weeks from B maximum, how- 
ever the spectral features of SN 2002ap remain broader than 
those of SN 2007gr This suggests that the differences between 
SNe 2007gr and 2002ap result from the differences in their den- 
sity profiles, i.e. a steeper density gradient for SN 2007gr and 
hence a smaller velocity range at which the lines can form. The 
photospheric velocity of SN 2007gr is also similar to that of 
SN 19941 about one week after B maximum, implying that the 
kinetic energy to unit mass should be also similar at this epoch. 

An interesting point to note is that the measured expan- 
sion velocities for Na I 5890, 5896 A about two weeks post- 
maximum show an increase for all Ic SNe in Fig{T5] with the 
exception of SN 2004aw, followed by a plateau. A possible ex- 
planation for this behaviour may be the contribution of another 
line. A likely possibility is the contrib ution from He 1 /t 5876 A 
which has been previously discussed bv lClocchiatti et akl (Il996l) . 
although as noted previously other strong lines of He are not ob- 
vious in the optical spectra. In addition to the possibility of con- 
tamination from He 1, the recombination of Na 11 to Na 1 may 
be an important mechanism which contributes to the strength of 



the Na I line in different velocity ranges. The observed behaviour 
may also be explained by the departure from LTE to NLTE, again 
affecting the strength of the Na ID feature. 

5.5. Carbon monoxide emission 

Ca rbon monoxide em i ssion in a supernova was first reported 
bv lCatchpole &Glassl ( |1987|) in SN 1987 A. Spectroscopic cov- 
erage of the first overtone (Av=2) sp anning days ~ 110-574 
was presented by iMeikle et al.l (1 19931) while coverage of the 
fundamental (Av =l) band emission spannin g d ays ~ 157-615 
was pr esented bv lBouchet & Danzigej(ll993h and lWooden et al.l 
(Il993h . Since then the first overtone of CO has been ob- 
served in numero us type II SNe with detections in the spec- 
tra of SN 199 5ad (Spvromilio & Leibundgu^[1996|), SN 1998S 
(iGerardv et al., 2000; Fassia et al. 20013 and SNe 1998dl & 
1999em (ISpvromiUo et al.ll200lb . 

With the sole exception of SN 1987 A, observations of the 
CO fundamental band (near ~4./im) have remained elusive, 
given the challenging nature of the observations. However, 
with space-based instrumentation, the CO fundamental has also 
been observed in several IIP SN e: SNe 2004dj, 2005af, 2004et 
(iKotak et al.ll2005l l2006l [20091) . Although CO emission in type 
II seems to be a common occu rrence, detections i n type Ic 
have been Hmited to SN 2000ew jGerardv et al.ll2002l) and now 
SN 2007gr 

The subject of molecular emission in SNe has been exten - 
sively reviewed by Petuchowski et al.'('l989^,'L epp et al. I (ll990l) . 
Dalgarno et al. (1998) and Gearhart et al. (19991), with particu- 
lar reference to SN 1987 A. The main formation mechanism is 
the radiative association of C and O atoms whereby photons are 
emitted and bound CO molecules are formed. With time, the flux 
of the CO emission will decrease as a result of destructive in- 
teractions with fast electrons which are produced by the y-rays 
from the decay of ^^Co. Additionally, the presence of ionised He 
that is microscopically mixed within the carbon-oxygen layers 
could destroy the CO molecules via charge transfer reactions: 

He+ + CO ^ C+ + O -H He (2) 

Thus molecular emission is a useful diagnostic for examining 
the conditions and degree of mixing in the ejecta of SNe. 

We see clear evidence for the presence of the first-overtone 
of CO emission from 2.29 fim to the red end of the spectrum 
at 2.4 /Ym. This feature first appears at H-70d post-B maximum 
which is the earliest reported detection of CO in any type of 
SN. In part, this might be due to the lack of suitable data for 
the majority of SNe. However, the early appearance of CO is in 
line with its relatively rapid evolution of SN 2007 gr compared to 
other SNe for which CO has been detected (e.g. SN 1987A) and 
is also consistent with the idea that the progenitor of SN 2007gr 
was a compact star. In earlier sections we have argued against 
the presence of large amounts of He in the ejecta of SN 2007gr, 
so the fact that CO is detected at all would seem to indirectly 
support this idea. Although one might invoke the possibilities 
that ionised He was n ot microscopically mixed with the un der- 
lying C/O rich layers dLepp et al.lll990l:lGerardv et al.ll2000l and 
references therein), or that a significant quantity of He was not 
ionised, both possibilities seem contrived in view of the obser- 
vation that strong spectral signatures due to He are not convinc- 
ingly detected over a time-span of more than a year. 

The CO feature grows in strength up to day H- 138 before fad- 
ing. In SN 1987A, a similar trend was observed, except that the 
intensity of the feature increased until ~ 200 d. This is believed 
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to be due to the non-uniform deposition of energy at early times 
as a result of the greater optical depth of y-rays in the iron core. 
After sufficient expansion of the ejecta by late times, the optical 
depth has decreased, and greater energy depositi on occurs in lay- 
ers ot her than the iron core, i.e. the C/O layer jLiu & Dalgarnd 
[1991) . Thus the energy deposition becomes progressively uni- 
form. 

Fig.[T2lshows that by day + 176, the CO emission has faded 
by a factor of ~ 1.5 compared to day H- 138 beyond 2.35 fim. This 



decrease in flux may represent the enhanced energy deposition 
rate into the C/O layers from fast electrons. 

As mentioned bv lGerardv et al.l (l2002h . the shape of the CO 
profile provides a lower limit to the temperature of the emitting 
region and so band profile modeling can be used to infer the ex- 
citation temperature. Since observations of CO emission in type 
Ic are rather limited, a comparison with the molecular emission 
in the spectra of SN 1987A is a good alternative. In Figure [16] 
we compare our /T-band spectrum of SN 2007gr at day + 137.7 
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to those of SN 1987A jMeikle et al.l 1 19891) and the type IIP 
SN 2004dj dKotak et al.ll2005h . The spectra have been scaled so 
that the initial slopes of the CO emissions near 2.24 yum are rela- 
tively well matched. In the + 137.7 day spectrum of SN 2007 gr, 
the CO feature decreases in flux at longer wavelengths in a sim- 
ilar manner to both SNe 1987 A and 2004dj. However, the first 
two band heads (labelled at -2.31 and 2.33 jum) of SNe 1987A 
and 2004dj can be clearly observed while that of SN 2007gr is 
not as distinct. LTE modelling of CO has shown that the two 
band heads have compa r able strengths at temperatures < 2000 K 
(ISharp & Hoflichlll989t lETu. Dalgarno & Lepplll992h . We thus 
conclude that the temperature of the CO emitting region in 
SN 2007gr is greater than 2000 K. 

An approximate velocity of the CO emission can be ob- 
tained by comparing our spectra of SN 2007 gr with the higher 
resolution spectrum of SN 1987 A. The slope of the 2-0 R- 
branch band head indicates that the expansion velocity in 
SN 2007 gr is approximately sim ilar to SN 1987 A (1800- 
2000 km s ' ISpvromilio et al.ll988|). which was a lso the case for 
SNe 1998dl and 1999em dSpvromilio et al.ll200ll) . 

CO was still detected in SN 1987A as late as 574 d after ex- 
plosion (iMeikle et al.ll 19931) . whereas it had disappeared by day 
355 in the type Iln SN 1998S dderardv et al . 2000). 

For SN 2007gr, we notice a fading of the feature by day 
+ 176, and we do not detect it in the + 377 d spectrum. However, 
we cannot entirely rule out the presence of the CO overtone 
since t he S/N is not sufficiently high. As suggested by lLiu et"aD 
( Il998l) . the disappearance of molecular emission may be due to 
the rapid decline in temperature as a result of the efficient cool- 
ing by CO at later times (decrease in heating rate and increasing 
departure from LTE). 

Aside from the diagnostic constraints e.g. on microscopic 
mixing, that are provided by CO, the presence of molecules in 
the ejecta provides an additional source of cooling. This may be 
conducive to the formation of dust in the ejecta. One might ex- 
pect the combination of molecules and low velocities observed 
in SN 2007gr to result in some condensation of dust grains, 
particularly graphite grains for the carbon-rich progenitor of 
SN 2007gr jValenti et al.ll2008al) . However, at the epoch of our 
last near-IR spectrum (~1 year), we do not see any clear evidence 
of a rising continuum in the /iT-band. The //-band light curve 
does not show any sign of flattening but appears to follow the 
^^Co decay rate reasonably well (see Fig.|4|i. The H - K colour 



evolution of SN 2007gr argues against dust formation with a rel- 
atively blue H-K colour of ~- 1 . 1 mag at ~ H- 400 d. Additionally, 
dust formation should result in a blueshift in the peaks of emis- 
sion lines which increases with time. We observe no such effect 
in the evolution of the emission features in the optical or near-IR 
spectra of SN 2007gr (see Fig.|2]&[l2l respectively). 

In relation to other core-collapse SNe types, IR sig- 
natures of dust condensation were found in the ejecta of 
the type IIP SN 1987A as early as 100 d post-explosion 



dSuntzeff & Bouc hd Il990h . Detections were also reported at 



e.g. feouchet & Danzigeill9"93h . Signatur es of dust were eviden t 
in flie type IIP SN 1999em after ~ 450 d dElmhamdi et alJl2003[) 
and also in SN 20 04et starting at ~300davs dSahu et al. 112006 : 
'Kotak et al ."20091). Additionally, dust formation in the type lb 
SN 199 01 was reported to have commenced as early as ~ 250 d 
(lEknhamdi et al.l2004l) . However, at the latest epochs under con- 
sideration here, we can conclusively state that signatures of dust 
formation in the ejecta of SN 2007gr are not apparent. Yet, given 
the rapid evolution of the SN, one might have expected signa- 
tures due to dust at epochs earlier than those for SN 1987 A, and 
other type II-P SNe. 

5.6. Estimate of the Fe mass 

The very late-time near-IR spectrum at H- 376.7 d (Fig. [12]) rep- 
resents the only one of its kind for type Ic SNe. 

As SN 2007gr evolves into its super-nebular phase, the rapid 
decrease in the electron temperature and density of the ejecta 
implies that strong [Fe II] and [Si I] lines in the J and H bands 
should be present. This is borne out by the observations with the 
prominent feature at ~1.65 fim either arising from [Fe II] 1.6435 
fim (a''F-a''D) or [Si I] 1.6454 fim, or even a blend of both. To 
further add to the complexity of this scenario, the intensity ratio 
between the lines of these ions (i.e. [Fe 11] 1.6454, 1.6068 fim 
and [Si I] 1.6435, 1.5994^/m) are expected to be similar at high 
densities (see ,01iva. Moorwood and Danziger 1987). 

However, if this feature is mainly due to [Fe II] then we 
should be able to detect also the 1.257 //m (a^D-a'*D) transi- 
tion, which shares the same upper level with the 1.6435 jum 
emission line. In our late-time spectrum, we do not clearly de- 
tect the 1.257 fim line which should, in fa ct, have an intensity 
-1.35 times greater than the 1 .6435 fim Une dMeikle etalJl989l) . 
Alternatively, if we assume [Si I] dominates the contribution to 
the 1 .65 fim feature, we expect to observe the strongest near-IR 
transition of [Si I] at 1.0991 fim. Thus the feature at ~ 1.1 fim 
may be mainly attributed to [Si I]. The identification of [Si I] is 
supported by the absence of a weaker [Fe II] emission at 1 .68 
fim (~ 25 % the flux of [Fe II] 1.6454yum). 

For the purpose of estimating an upper limit to the mass of 
^^Fe ejected in SN 2007 gr , we have used the earlier spectrum at 
+ 175 d where the identification of the 1.65 fim feature with [Fe 
II] is more secure and the S/N is much higher 

The quantity of iron in the ejecta can be derived if we know 
the tem perature, density and emissivity of the [Fe II] emission 
regions dGraham et al.lll98^ . Therefore models which consider 
the ionisation structure of the SN ejecta are necess ary to deter- 
mine an accurate iron mass at this epoch. Graham et al.l dl986h 
related the 1 .6435 micron line flux to the iron mass (Mfe) for the 
Type lb SN 1983N according to the equation; 



350 d (Meikle et al.ll 19931) and ~500 days after the outburst (see 



FiMAiim = 77-75 — —exp{-TH) (3) 



npeii npe mfe 
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where A denotes the Einstein spontaneous decay probability 
(0.00569 s"'; Nussbaumer & Storey 1982), /iv is the energy of 
the transition (i.e. a4Fc)/2-a4D7/2), and d is the distance to the 
SN. The optical depth of the lin e is represented by th , and mpe 
denotes the mass of an iron atom. lGraham et an(ll986l) estimated 
that the ejecta of SN 1983N contained -0.3 M© of iron. 

To obtain an estimate for SN 2007gr, we have adopted sim- 
ilar physical parameters. Assuming LTE conditions, the forbid- 
den lines are excited at temperatures ~6500K, and so na4D^^Jnfe 
~ 0.03. This assumption is valid for our 176 d spectrum since 
electron densities ar e typically > 10* cm"^ before ~ 200 days 
(lGrahametal.1119861) . 

Since our spectrum was observed ~ half a year from B max- 
imu m, the ionisation fraction npeii /nfe can be approximated as 
0.1 (iMeverotll 1 9801: lAxelrodl 1 980l) . Taking into account the dis- 
tance to SN 2007gr and a negligible extinction for the 1 .6435 jum 
line in the H band (see Table [U, the integrated line flux of ~6 
xlO"'*" erg cm^^s"' implies that 0.14 + 0.08 Mq of Fe was syn- 
thesised in SN 2007 gr Since the quantity of ^^Ni derived from 
the bolometric light curve is 0.076 + 0.010 Mq, the iron mass is 
a plausible upper limit. The error in the estimate of Fe accounts 
for the uncertainties in the distance to SN 2007gr, the extinction 
and also on the atomic data. 

6. Summary 

We have presented an extensive set of optical and near-IR obser- 
vations of the type Ic SN 2007gr. In order to put SN 2007gr into 
context, we compared our data with a sample of relatively well- 
monitored type Ic SNe that span a range in observed properties. 
The light curves of SN 2007gr best resemble those of the broad- 
lined ('hypernova') type Ic SN 2002ap at early times, peaking 
with a B band magnitude of Mb - -16.8. However, the peak 
absolute optical magnitudes lie towards the fainter end of the Ic 
SNe that we considered here. 

The narrow spectral features of SN 2007gr compared to 
typical SNe Ib/c, are indicative of low expansion velocities, 
and facilitate line-identification which is usually problematic in 
type Ic SNe due to severe line-blending. The similarity with 
SN 2002ap does not extend to the spectroscopic observations, 
with SN 19941 generally providing a better match. 

Our nebular phase spectra reveal a systematic blue-shift in 
the profiles of forbidden [O I] and Mg I], while those of the 
permitted lines of the same elements are centred at their rest 
wavelengths. This implies that the forbidden and permitted lines 
appear in different regions of the ejecta, as expected. The sys- 
tematic blue-shift of the forbidden lines is mostly likely due to 
the absorption of photons emitted on the far side of the SN by 
the intervening, higher opacity material. 

We found low Mg I]/[0 I] ratios for SN 2007gr during the 
early nebular phase supporting the previous suggestion that the 
progenitor was a carbon-rich star ( Valenti et al. 2008a). This 
finding, taken together with the comparisons with SN 2002ap, 
and the relatively rapid spectral evolution e.g. the early appear- 
ance of CO, suggest a highly-evolved, compact progenitor. From 
modelling of the quasi-bolometric light curve, we estimate that 
a modest 2-3.5 of material was ejected by SN 2007 gr at a 
kinetic energy of 1-4 foe, and 0.076 + 0.020 Mq of ^''Ni was 
synthesised in the explosion - with the latter value being very 
similar to that found for SN 2002ap. 

We clearly detected the first-overtone of the CO molecule 
and were able to track its evolution. It first appeared about day 7 1 
- the earliest recorded appearance of this molecule and persisted 
for at least 100 d. 



The proximity of SN 2007 gr, one of the nearest SNe Ic in 
recent years, has provided us with an ideal target for intensive 
follow-up. Indeed, SN 2007 gr is the first type Ic SN for which 
near-IR spectra have been obtained at all phases. We suggest that 
SN 2007gr is an ideal template for future SNe Ib/c studies. 
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Table A.l. U, B, V, R and / magnitudes of the sequence stars in the field of NGC 1058. The errors reported in parentheses are the 
rms of the measurements of the different photometric nights. 



Star 



U 



1 18.531 (0.020) 

2 15.039 (0.010) 

3 16.904 (0.020) 

4 17.242 (0.020) 

5 18.158 (0.030) 

6 17.920 (0.050) 

7 16.492 (0.030) 

8 15.504 (0.012) 

9 14.186(0.010) 

10 13.916(0.010) 



B 

17.232 (0.009) 
15.052 (0.005) 
16.400 (0.011) 
17.230 (0.009) 
17.684 (0.014) 
16.913 (0.018) 
15.544 (0.007) 
15.431 (0.007) 
14.247 (0.010) 
13.876 (0.006) 



V 

16.042 (0.010) 
14.514 (0.010) 
15.544 (0.013) 
16.578 (0.010) 
16.857 (0.014) 
15.804 (0.020) 
14.431 (0.007) 
14.799 (0.008) 
13.675 (0.005) 
13.234 (0.010) 



R 

15.277 (0.010) 
14.184 (0.010) 
15.047 (0.011) 
16.200 (0.008) 
16.384 (0.015) 
15.111 (0.013) 
13.835 (0.005) 
15.125 (0.010) 
13.339 (0.009) 
12.858 (0.010) 



__/__ 

14.655 (0.015) 
13.876 (0.012) 
14.614 (0.009) 
15.832 (0.010) 
15.961 (0.008) 
14.476 (0.012) 
13.300 (0.010) 
14.770 (0.008) 
13.010(0.013) 
12.457 (0.011) 
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Table A.2. Optical Photometry of SN 2007gr " 
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Table A.2. continued. 
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4 


19/10/07 


54392.60 


+55.6 










15.984 


+ 


0.028 


15.162 ± 


0.025 


14.680 ± 


0.022 


14.060 + 


0.020 


1 


21/10/07 


54394.52 


+57.5 


16 


571 ± 





038 


16.058 


+ 


0.032 


15.201 ± 


0.026 


14.743 ± 


0.023 


14.127 ± 


0.025 


4 


23/10/07 


54396.52 


+59.5 










16.073 


+ 


0.023 


15.235 ± 


0.024 


14.772 ± 


0.021 


14.133 ± 


0.022 


4 


24/10/07 


54398 48 


+61.5 


16 


591 ±0 


037 


16.142 


+ 


0.029 


15.267 ± 


0.022 


14.853 ± 


0.024 


14.168 ± 


0.025 


4 


29/10/07 


54402.58 


+65.6 


16 


690 ± 





043 


16.168 


+ 


0.031 


15.319 ± 


0.025 


14.906 ± 


0.024 


14.216 ± 


0.023 


4 


30/10/07 


54403.58 


+66.6 










16.175 


± 0.029 


15.352 ± 


0.021 


14.905 ± 


0.020 


14.238 ± 


0.020 


5 


02/11/07 


54406.70 


+69.7 


16 


662 ± 





042 


16.209 


+ 


0.026 


15.379 ± 


0.024 


14.973 ± 


0.025 


14.241 ± 


0.021 


4 


06/11/07 


54410.70 


+73.7 


16 


/4U ± 





036 


16.295 


+ 


0.025 


15.458 ± 


0.021 


15.046 ± 


0.023 


14.351 ± 


0.023 


4 


13/11/07 


54417.58 


+80.6 










16.295 


+ 


0.031 


15.600 ± 


0.022 


15.133 ± 


0.023 


14.518 ± 


0.025 


1 


14/1 1/07 


54418.56 


+81.6 


16 


O 1 1 1 





051 


16.378 


+ 


0.041 


15.579 ± 


0.030 


15.177 ± 


0.027 


14.480 ± 


0.032 


4 


21/11/07 


54425.63 


+88.6 










16.472 


+ 


0.020 






15.215 ± 


0.020 






5 


27/11/07 


54432.44 


+95.4 


17 


027 ± 





045 


16.581 


+ 


0.037 


15.824 ± 


0.022 


15.388 ± 


0.026 


14.724 ± 


0.028 


4 


04/12/07 


54439.51 


+ 102.5 










16.607 


+ 


0.040 


15.923 ± 


0.031 


15.427 ± 


0.025 


14.825 ± 


0.023 


1 


05/12/07 


54440.45 


+ 103.5 










16.665 


+ 


0.035 


15.938 ± 


0.025 


15.442 ± 


0.023 


14.922 ± 


0.024 


1 


06/12/07 


54440.53 


+ 103.5 


17 


094 ± 





040 


16.657 


+ 


0.032 


15.935 ± 


0.023 


15.506 ± 


0.021 


14.859 ± 


0.023 


4 


10/12/07 


54444.51 


+ 107.5 










16.745 


+ 


0.028 


16.039 ± 


0.026 


15.519 ± 


0.021 


15.000 ± 


0.021 


5 


10/12/07 


54444.54 


+ 107.5 


17 


203 ± 





060 


16.761 


± 0.032 


16.028 ± 


0.022 


15.570 ± 


0.022 


15.004 ± 


0.025 


4 


16/12/07 


54450.62 


+ 113.6 


17 


188 ± 





070 


16.830 


+ 


0.040 


16.122 ± 


0.030 


15.626 ± 


0.027 


15.058 ± 


0.028 


2 


25/12/07 


54460.49 


+ 123.5 










16.970 


+ 


0.023 


16.326 ± 


0.023 


15.709 ± 


0.026 


15.206 ± 


0.026 


5 


28/12/07 


54463.31 


+ 126.3 










17.030 


+ 


0.029 






15.781 ± 


0.026 


15.278 ± 


0.026 


5 


28/12/07 


54463.39 


+ 126.4 


17 


327 ± 





060 


17.022 


+ 


0.047 


16.346 ± 


0.025 


15.781 ± 


0.027 


15.226 ± 


0.024 


1 


03/01/08 


54468.50 


+ 131.5 










17.077 


+ 


0.031 


16.476 ± 


0.024 


15.803 ± 


0.021 


15.343 ± 


0.021 


5 


06/01/08 


54472.40 


+ 135.4 










17.180 


+ 


0.045 


16.495 ± 


0.031 


15.935 ± 


0.025 


15.442 ± 


0.025 


1 


09/01/08 


54474.51 


+ 137.5 










17.238 


+ 


0.037 


16.568 ± 


0.024 


15.922 ± 


0.022 


15.419 ± 


0.020 


5 


13/01/08 


54479.40 


+ 142.4 


17 


808 ± 





055 


17.239 


+ 


0.040 


16.725 ± 


0.037 


15.959 ± 


0.030 


15.538 + 


0.027 


3 


26/01/08 


54492.35 


+ 155.4 


17 


871 ± 





131 


17.392 


+ 


0.046 


16.853 ± 


0.028 


16.150 ± 


0.025 


15.810 + 


0.026 


8 


28/01/08 


54494.36 


+ 157.4 










17.477 


+ 


0.047 


16.967 ± 


0.024 


16.161 ± 


0.020 


15.829 ± 


0.020 


5 


10/02/08 


54507.45 


+ 170.5 










17.671 


+ 


0.040 






16.371 ± 


0.027 


16.168 ± 


0.025 


5 


29/08/08 


54707.97 


+371.0 
















20.300 ± 


0.200 


19.770 ± 


0.200 


19.620 ± 


0.200 


9 


02/09/08 


54711.62 


+374.6 
















20.400 ± 


0.500 


19.800 ± 


0.500 


19.540 ± 


0.400 


8 


30/09/08 


54740.49 


+403.5 
















20.696 ± 


0.500 


20.205 ± 


0.400 


19.861 ± 


0.600 


8 



a The S-corrected optical photometry of SN 2007gr The reported errors (in parentheses) are a result of the uncertainty due to the PSF fitting of 
the SN magnitude and to the background contamination (computed using the artificial-star experiment). 
b Relative to the epoch of B maximum JD = 2454337.0 ± 1.0. 



1= Ekar: 1.82 m Copernico Telescope, INAF - Osservatorio di Asiago, Mt. Ekar, Asiago (Italy) + AFOSC. 

2= TNG: 3.5 m Telescopio Nazionale Galileo, La Palma, Canary Isl. (Spain) + DOLORES. 

3= NOT: 2.5 m Nordic Optical Telescope, La Palma, Canary Isl. (Spain) + ALFOSC. 

4= LT: 2.0 m Liverpool Telescope, La Palma, Canary Isl. (Spain) + RATCAM. 

5= Wendelstein: 0.8 m Telescope at the Wendelstein Observatory, (Munich, Germany) + MONICA. 

6= S70: 0.70 m Sternberg Astronomical Institute Telescope, Moscow (Russia). 

7= G. D. Cassini: 1.52 m Telescope at the Loiano Observatory (BOL., Italy) + BFOSC. 

8= CA: 2.2 m Telescope at the Calar Alto Observatory (Spain) + CAFOS. 

9= Gemini North: 8. 1 m Telescope, Mauna Kea (Hawaii) + GMOS-N. 
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Table A.3. Near-Infrared Photometry of SN 2007gr " 



U.T. Date 


J.D. 
2,400,000+ 


Phase * 
(days) 




J 






H 






K 






Instrument 


22/08/07 


54335.74 


-1.3 


12 


570 ± 





100 














1 


02/09/07 


54347.50 


+ 10.5 


12 


471 ± 





030 


12.306 ± 


0.030 


12 


077 ± 





030 


2 


05/09/07 


54349.50 


+ 12.5 


12 


518± 





020 


12.307 ± 


0.020 


12 


154 ± 





030 


3 


18/09/07 


54361.62 


+24.6 










12.760 ± 


0.040 










4 


19/09/07 


54363.60 


+26.6 


13 


069 ± 





020 


12.723 ± 


0.030 


12 


651 ± 





030 


5 


08/10/07 


54382.60 


+45.6 


13 


733 ± 





040 


13.434 ± 


0.040 


13 


377 ± 





040 


1 


10/10/07 


54383.39 


+46.4 










13.440 ± 


0.050 










4 


14/10/07 


54388.60 


+51.6 


13 


913 ± 





040 


13.507 ± 


0.040 


13 


681 ± 





040 


1 


17/10/07 


54391.60 


+54.6 


14 


020 ± 





040 


13.657 ± 


0.040 


13 


674 ± 





040 


1 


24/10/07 


54398.60 


+61.6 


14 


239 ± 





040 


13.807 ± 


0.040 


13 


877 ± 





040 


1 


28/10/07 


54402.60 


+65.6 


14 


411 ± 





040 


13.911 ± 


0.040 


13 


957 ± 





040 


1 


03/11/07 


54407.32 


+70.3 










14.030 + 


0.030 










4 


20/11/07 


54424.44 


+87.4 










14.490 ± 


0.030 










4 


05/12/07 


54439.70 


+ 102.7 


15 


396 ± 





040 


14.795 ± 


0.030 


15 


309 ± 





030 


3 


16/12/07 


54450.70 


+ 113.7 


15 


860 ± 





040 


15.120 ± 


0.030 


15 


415 ± 





030 


3 


20/12/07 


54454.35 


+ 117.4 










15.240 ± 


0.040 










4 


09/01/08 


54474.20 


+ 137.2 










15.590 + 


0.040 










4 


09/01/08 


54474.72 


+ 137.7 


16 


306 ± 





100 


15.552 + 


0.100 


15 


603 ± 





100 


5 


16/02/08 


54512.22 


+ 175.2 










16.250 ± 


0.040 










4 


20/02/08 


54516.23 


+ 179.2 










16.320 ± 


0.040 










4 


05/09/08 


54715.13 


+378.1 


19 


214 ± 





060 














4 


06/09/08 


54716.14 


+379.1 










18.196 + 


0.029 










4 


08/09/08 


54718.04 


+381.0 














19 


280 ± 





200 


4 



a The reported errors (in parentheses) are a result of the uncertainty due to the PSF fitting of the SN magnitude and to the background 
contamination (computed using the artificial-star experiment). 
b Relative to the epoch of B maximum (JD = 2454337.0 ± 1.0). 

1= CI: 1.1 m Telescope AZT-24, Campo Imperatore (Italy) + SWIRCAM. 

2= NOT: 2.5 m Nordic Optical Telescope, La Palma, Canary Isl. (Spain) + NOTCAM. 

3= TNG: 3.5 m Telescopio Nazionale Galileo, La Palma, Canary Isl. (Spain) + NICS. 

4= Gemini North: 8.1 m Telescope, Mauna Kea (Hawaii) + NIRI. 

5= UKIRT: 3.8 m Infrared Telescope, Mauna Kea (Hawaii) + WFCAM. 



Hunter e? al: SN 2007gr 
Table A.4. Log of Optical Spectroscopy for SN 2007gr. 



U.T. Date 


JD 


Epoch " 


Instrumental 


Resolution 


Exptime 


Range 


Reference ' 




(2,400,000+) 


(days) 


configuration 


(A)* 


(s) 


(A) 




2007 Aug 17 


54329.7 


-7.3 


TNG+DOLORES+gm.LRS+LRB 


12 


900 


3200-7710 


1 


2007 Aug 17 


54329.7 


-7.3 


NOT+ALFOSC+gm.4 


12 


250 


3400-9060 


2 


2007 Aug 18 


54330.5 


-6.5 


Ekar+AFOSC+gm.4,gm.2 


15-25 


1200 


3490-9980 


1 


2007 Aug 20 


54332.7 


-4.3 


NOT+ALFOSC+gm.4 


12 


900 


3240-9110 


1 


2007 Aug 20 


54332.7 


-4.3 


TNG+DOLORES+gm.LRS+LRB 


12 


200 


3180-9400 


1 


2007 Aug 21 


54333.7 


-3.3 


NOT+ALFOSC+gm.4 


12 


300 


3170-9100 


1 


2007 Aug 22 


54334.7 


-2.3 


NOT+ALFOSC+gm.4 


12 


200 


3190-9110 


1 


2007 Aug 23 


54335.7 


-1.3 


NOT+ALFOSC+gm.4 


12 


300 


3130-9100 


1 


2007 Aug 26 


54338.5 


+1.5 


Pennar +B&C+300tr/mm 


25 


600 


3270-7810 


1 


2007 Aug 26 


54339.7 


+2.7 


NOT+ALFOSC+gm.4 


12 


600 


3190-9100 




2007 Aug 31 


54344.5 


+7.5 


Pennar +B&C+300tr/mm 


25 


1200 


3450-7800 


1 


2007 Sep 07 


54350.7 


+ 13.7 


NOT+ALFOSC+gm.4 


12 


600 


3200-8940 


1 


2007 Sep 10 


54353.6 


+16.6 


NOT+ALFOSC+gm.4 


12 


600 


3200-9060 


1 


2007 Sep 14 


54357.6 


+20.6 


Loiano+BIFOSC+gm.4,5 


15-25 


900 


3500-9780 


1 


2007 Sep 14 


54357.5 


+20.5 


Ekar+AFOSC+gm.4,2 


15 


1200 


3550-9830 


2 


2007 Sep 15 


54358.6 


+21.6 


NOT+ALFOSC+gm.4 


12 


900 


3070-9270 


2 


2007 Sep 21 


54364.6 


+27.6 


Pennar +B&C+300tr/mm 


25 


3000 


3270-7830 


2 


2007 Sep 22 


54365.7 


+28.7 


NOT+ALFOSC+gm.4 


12 


900 


3180-9070 


2 


2007 Sep 25 


54368.6 


+31.6 


NOT+ALFOSC+gm.4 


12 


900 


3200-9100 


2 


2007 Sep 30 


54373.7 


+36.7 


NOT+ALFOSC+gm.4 


12 


900 


3240-9070 


1 


2007 Oct 09 


54382.6 


+45.6 


Ekar+AFOSC+gm.4,2 


15-25 


2400 


3490-10050 


2 


2007 Oct 19 


54392.6 


+55.6 


Ekar+AFOSC+gm.4,2 


15-25 


1800 


3540-9860 


2 


2007 Nov 13 


54417.5 


+80.5 


Ekar+AFOSC+gm.4,2 


15-25 


2700 


3490-9980 


2 


2007 Dec 05 


54440.4 


+103.4 


Ekar+AFOSC+gm.4,2 


15-25 


2700 


3490-10070 


2 


2007 Dec 15 


54450.5 


+113.5 


TNG+DOLORES+VHR+LRS 


5 


1800 


4660-6640 


2 


2008 Jan 06 


54472.3 


+135.3 


Ekar+AFOSC+gm.4,2 


15-25 


1800 


3490-10180 


2 


2008 Jan 29 


54495.4 


+158.4 


CA+CAFOS+gm.b200,r200 


12 


2400x2 


3210-10630 


2,3 


2008 Sep 02 


54711.5 


+374.5 


CA+CAFOS+gm.g200 


12 


2700x3 


3670-10460 


2 


2008 Sep 03 


54712.5 


+375.5 


GN+GMOS+R150 


12 


5400 


4000-11000 


2 



a Relative to the epoch of B maximum (JD = 2,454,337.0 ± 1.0). 

b Derived from F WHM of n ight-sky emission lines. 

c References: ni lValenti e t al. (2008a) [2] This paper r31 IValenti et"al] ( 120091) . 



TNG = 3.5 m Telescopio Nazionale Galileo, La Palma, Canary Isl. (Spain) + DOLORES. 

Ekar = 1.82 m Copernico Telescope, INAF -Osservatorio di Asiago, Mt Ekar, Asiago (Italy) + AFOSC. 

NOT = 2.5 m Nordic Optical Telescope, La Palma, Canary Isl. (Spain) + ALFOSC. 

Pennar = 1.22 m Telescope, Asiago Pennar, (Italy) + Boiler & Chivens (B&C) spectrograph. 

Loiano = 1.52. m Telescope, INAF -Bologna Observatory, Bologna (Italy) + BIFOSC. 

CA = 2.2 m Telescope at Calar Alto Observatory (Spain) + CAFOS. 

GN = 8. m Telescope at Gemini North, Hawaii + GMOS-N. 
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Table A.5. Log of Near-Infrared Spectroscopy of SN 2007gr. 



U.T. Date 


J.D 


Epoch " 


Standard star 


t,,p/band SN " 


Instrumental 


Resolution 




(2,400,000+) 


(days) 


(Spectral type) 


(s) 


configuration ^ 


(A)"* 


2007 Sep 09 


54352.5 


+ 15.5 


HIP 1603 (AOVs) 


300 


TNG+NICS+gm.IJ,HK 


15-30 


2007 Sep 18 


54362.1 


+25.1 


HD 19301 (F3V) 


300 


GN+NIRI+JHK 


20-30 


2007 Oct 05 


54378.7 


+41.7 


HIP 1603 (AOVs) 


300 


TNG+NICS+gm.IJ,HK 


15-30 


2007 Oct 10 


54383.9 


+46.9 


HIP 12181 (F5V) 


300 


GN+NIRI+JHK 


20-30 


2007 Nov 03 


54407.8 


+70.8 


HD 17119 (F5V) 


300 


GN+NIRI+JHK 


20-30 


2007 Nov 20 


54424.9 


+87.9 


HD 19301 (F5V) 


300 


GN+NIRI+JHK 


20-30 


2007 Dec 16 


54451.5 


+ 114.5 


HIP 10559 (AOV) 


600 


TNG+NICS+gm.IJ,HK 


15-30 


2007 Dec 20 


54454.8 


+ 117.8 


HD 17119 (F5V) 


840 * 


GN+NIRI+JHK 


20-30 


2008 Jan 09 


54474.7 


+ 137.7 


HD 17119 (F5V) 


840 


GN+NIRI+JHK 


20-30 


2008 Feb 16 


54512.7 


+ 175.7 


HD 17119 (F5V) 


1120 


GN+NIRI+JHK 


20-30 


2008 Sep 04 


54713.6 


+376.6 


HIP 12719 (B3V) 


1800 


GN+NIRI+JHK 


20-30 



a Epoch with respect to the estimated i?-band maximum (ID 2,454,337.0 ± 1.0). 
b The quantity t„,,/band is the on-source integration time for each band. 

* 1050 s exposure for the / band. 
c TNG= 3.5 m Telescopio Nazionale Galileo, La Palma, Canary Isl. (Spain), GN= 8.1m Telescope at Gemini North, Hawaii. 
d The instrumental resolution for the Gemini spectra was derived as the range of Gaussian FWHMs measured for a number of unblended arc 
lines. The resolution of the TNG spectra have been estimated using the FWHM of the night sky lines. 



